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A model of alloy is considered which includes both quenched disorder in the electron 
subsystem ("alloy" subsystem) and electron-phonon interaction. For given approximate 
solution for the alloy part of the problem, which is assumed to be conserving in Baym's 
sense, we construct the generating functional and derive the Eliashberg-type equations 
which are valid to the lowest order in the adiabatic parameter. The renormalization of 
bare electron-phonon interaction vertices by disorder is taken into account consistently 
with the approximation for the alloy self-energy. For the case of exact configurational 
averaging the same set of equations is established within the usual T-matrix approach. 
We demonstrate that for any conserving approximation for the alloy part of the self- 
energy the Anderson's theorem holds in the case of isotropic singlet pairing provided 
disorder renormalizations of the electron-phonon interaction vertices are neglected. Tak- 
ing account of the disorder renormalization of the electron-phonon interaction we analyze 
general equations qualitatively and present the expressions for Tc for the case of weak 
and intermediate electron-phonon coupling. Disorder renormalizations of the logarithmic 
corrections to the effective coupling, which arise when the effective interaction kernel for 
the Cooper channel has the second energy scale, as well as the renormalization of the 
dilute paramagnetic impurity suppression are discussed. 

Keywords: Electron-Phonon Interaction, Quenched Disorder, Superconductivity, Gener- 
ating Functional Approach 

1. Introduction 

A problem of how disorder influences a supercondjjntive transition temperature Tc, 



first ^tudied in pioneering papers by Anderson JJr and_by Abrikosov and Gor'- 
kov, 13 draws a great deal of physicists attention.oBlj IiliJ'B'EiJ However, even if a usual 



mechanism of the superconductivity owing to electron-phonon interaction (EPI) 
and isotropic s-type pairing are concerned — it is this subject present paper deals 
with — one cannot regard the problem as being solved completely. The situation 
where disorder scattering is weak or the concentration of impurities (both magnetic 
and nonmagnetic) is small was studied rigQrously.&Q'Q 

As for strong disordering, in early work this problem was considered within 
the framework of semiphenomenological approaches which were called for explaining 
qualitatively empirical rules (e.g., Miedema rulesEj'cjEj) determining supercondiic- 
tive transition temperature in alloys. In a pioneering paper by Anderson et alJlj 
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a mechanism of universal Tc-degradation was proposed which was based on the 
enhancement of screened Coulomb pseudopotential due to disordering effects in 
compounds with A-15 structure. Further a great number of works appeared touch- 
ing interrelations between the phonon-induced superconducti vity an d Anderson's 
localization, and other effects of strong disorder, see papersJl3'OE3iijc3 and also a 
recent review by Belitz and Kirkpatricla where mapping on the non-linear matrix 
(T-model was used when treating the Anderson's localization and related problems. 

Recently a new impetus has appeared for studying the influence of disorder on 
Tc. This impetus comes mainly from discussing mechanisms of superconductiv- 
ity in novel Cu-oxide based compounds (see, e.g., papers by Abrikosov,S where 
the influence of non-magnetic impurities on Tc is discussed in the case of strongly 
quasimomentum-dependent singlet order parameter). These substances are sup- 
posed to be strongly correlated systems which may be described in the simplest 
way by an "aUpy analogy" within some approaches, as for instance Hubbard-IILaxt 
proximation E3 or static approximation within functional integral approachEi'EafiaO 
where this analogy emerges in a rather natural way. In particular, the renormaliza- 
tion of the phonon-induced Tc was discussed in Ref.c2 where the use was made of 
Hubbard-III approximation when treating the effects of strong electron on-site re- 
pulsion, however only a part of contributions was taken into account. Although the 
approximations of Hubbard-III type are considered now as not being quite adequate, 
from the formal point of view they have many similarities with mode r n anproa . c hes 
to highly correlated systems, which are based on the d ^ oo limit iESEjOHE^ the 
electron self-energy possesses strong energy dependence in the normal phase, while 
the quasimomentum dependence is absent. 

Thus, long time after creating BCS theory the classical problem of the disorder 
influence on Tc attracts attention of researchers. 

At the same time, existence of large number of works and approximations has 
lead to that it is difficult to compare results of various approaches between each 
other. Even the order parameter is understand in different ways. Therefore it is 
instructive to study the problem in a general form by which we mean the follow- 
ing. Given an approximate self-energy expression for the averaged electron Green's 
function in the normal phase and without electron-phonon interaction, we search 
for the equations determining the order parameter in the superconductive phase for 
the electron-phonon pairing mechanism, which are consistent (in some still unde- 
termined sense) with the approximation. 

In the present paper we study systematically the question of how strong quenched 
substitutional disorder (components of an alloy are assumed to be nonmagnetic) in- 
fluences Tc, leaving aside an interesting in its own right and fascinating problem of 
an interplay between superconductivity and localization phenomenon (the latter is 



known may occure in disordered systemsEj). We shall demonstrate that for a given 
approximation describing normal state of an alloy in the absence of EPI and satisfy- 
ing natural requirement of being thermodynamically consistent in Baym-Kadanoff 
sensa^O (later on rather the terms "0-derivable" or "conserving" will be used in- 
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terchangeably when referring to such approximations), and with suitable meaning 
of the term "consistent" , the problem can be solved in an unique way. We propose 
a procedure which enables constructing the Eliashberg-type equations, which are 
valid to the lowest order in the adiabatic parameter and properly take into account 
the disorder renormalization of EPI vertices. The structure of these equations is 
the same for all ^-derivable approximations and such approximations, as Virtual 
Crystal, Diffusional, Averaged T-Matrix and Coherent Potential approximations, 
etc., can be discussed within an unified scheme. However, the linearized version 
of the equations is so complicated that further simplifications are required and we 
establish also what kind of simplifications leads to the results previously obtained 
within some (in fact, within the most of already mentioned) approximations widely 
accepted in treating the "alloy" problems. 

We_^s]iQidd note that the problem in a similar form has been considered by 
BelitiaOH with the use of exact eigenfunction representation. However, in this 
approach the very definition of the order parameter differs from the traditional 
one, which exploits quasimomentum or coordinate representation, and establishing 
relations between these approaches is, in general, a very difficult problem. The 
approach being proposed in this paper seems to be more directly related with ob- 
servable characteristics, example of which is the quasimomentum dependence of the 
superconductive gap. In this connection, the traditional and often not specially 
discussed supposition, that the order parameter is constant over the Fermi sur- 
face in dirty metals and, therefore, the type of averaging, procedure is not crucially 
important, docs not, generally speaking, hold (see RefH3). 

The paper is organized as follows. In Section 2 we discuss the model. The gener- 
ating functional is constructed and the equations of the strong coupling theory are 
then derived. An alternative approach with the use of standard T-matrix formu- 
lation for the configurational averaging is given in Section 3 where we discuss also 
the relationships between the two approaches. In Section 4 we demonstrate that 
for the case of the isotropic singlet pairing owing to electron-phonon interaction the 
Anderson's theorem may be valid when formulated analogously to that for the BCS 
model, and we explicitly specify the conditions for the theorem to hold (all ear- 
lier treatments appear to be incomplete). In Section 5 reduced isotropic equations 
are derived in which the anomalous self-energy contribution coming from disorder 
scattering ( "disorder" contribution to the anomalous part of the self-energy) is elim- 
inated exactly. Section 6 deals with qualitative analysis of the reduced equations 
within certain approximations for the isotropic contributions of electron-phonon 
vertices and general expressions for Tc are presented there for the case of weak and 
intermediate electron-phonon coupling. We establish relations between effective 
electron-phonon coupling and generating function for arbitrary thermodynamically 
consistent single-site approximation for the alloy self-energy and write down ex- 
plicit expressions for renormalized electron-phonon coupling within a number of 
widely known single-site approximations in Section 7. In Appendix A we give a 
derivation of Ward-type identities. Appendix B deals with renormalizations of the 
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Tc-suppression by dilute paramagnetic impurities and in Appendix C some use- 
ful relations between generating functions for "disorder" self-energy and for vertex 
corrections are established within arbitrary single-site approximation. 

2. The Model, Generating Functional and Eliashberg-type Equations 

To begin with, we consider the model of alloy with substitutional disorder described 
by the Frolich-type model Hamiltonian 

H = He + He-ph + Hph, (2.1) 

where He is Hamiltonian of an electron subsystem, which in the Wannier represen- 
tation has the form 

^e = ^%ctcja, (2.2) 

hij = tij - fiSij + + EiSij, (2.3) 

and ctj. (da) is electron creation (annihilation) operator in the Wannier state with 
spin cr on site i; iy, non-random part of transfer integral; ^^j, random part of 
transfer integral and < >= (< ... > means configurational average with 
certain probability distribution); Si, on-site electron energy; Sij, usual Kronecker 
delta symbol and ^, the chemical potential. 

In the case of substitutional disorder of a general type both and Si are random 
variables which depend on kinds of atoms placed at the corresponding lattice sites 
and particular cases of alloy disorder may be obtained by appropriate choosing the 
quantities and Ei to be random or not. The so called diagonal disorder emerges 
when one sets all the to be identically zero and quantities Ei are random and 
depend on the sort of atoms at the corresponding lattice sites. If, on the contrary, 
are random, the Hamiltonian ( |2.2| ) describes the case of off-diagonal alloy disorder .Ej 

To complete the description of a disorder type one should fix distributions for the 



random quantities entering the Hamiltonian (2.2). For the case of random binary 



alloy model with diagonal disorder Eq. (2^) takes the form 



hij = Uj - fiS^j + eaSij + Vr]i6^j, (2.4) 
Prob(7yi = 1) = c, Prob(r/i ^ 0) = 1 - c, (2.5) 

with c being the concentration of B-typc component of the alloy in the host which 
is assumed to be of A- type, and V = Eb — ea is the scattering potential of B-type 
atoms measured relatively to the host potential. We will use this specific model of 
alloy in Sections 6 and 7 and in Appendices B and C. In the remaining part of the 
paper a discussion is rather general and does not rely upon any specific alloy model. 
The Hamiltonian of the phonon subsystem is of the following form 
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with 6q ^ (foq,s) being phonon creation (annihilation) operator with quasimomentum 
q and the branch index s, ti^q^s being the phonon frequencies and we assume that 
wq,s are non-random functions despite they are renormalized by all the interactions 
in our model. 

Finally, the Hamiltonian of the electron-phonon interaction reads 

J/e-p. = E E ^434'^.-^ E exp(-iqRO {h\ s + fo-q,.) , (2.7) 



ij.a l,s ' q 

where M^^^^ is matrix element of the electron-phonon interaction (EPI) and supposed 
to be independent of disorderJil 

Is) 

Deriving explicite expressions for M)- - requires some care since because the 
scattering by impurities is ellastic only in the coordinate frame moving with crystall 
lattice as was pointed out by Blount. l3 We will not write down the expressions 

is) 

for corresponding detalies and the results of calculations are given in 

Refs 



As is knownii one should treat cjq s as completely renormalized by all the rel- 
evant interactions when considering the phenomenon of superconductivity in such 
model, and effects of EPI on electron spectrum should be taken into account in the 
lowest (in fact, second) order in Mj*:^j for each disorder configuration. This actually 
corresponds to the lowest order of the pertubation theory in addiabatic parameter 
as was demonstrated by MigdaLcj Of course the two simplifying approximations 

is) 

about non-random character of M^.^ - and cjjq.s cannot be proven regorously if the 
real situation is kept in mind, and therefore they have a model character within the 
framework of the Frolich-type model. However, as we will see in what follows, even 
with these approximations the model remains complicated and interesting enough 
for studying. 

In the Gorkov-Nambu representationao (2^) and (2/7) take the form 



i/e^E'^'^cJraC,, (2.8) 
He^vh = T.T. ^43cJ^3C,^ E exp(-iqRO {bl^ + b-d.s) , (2.9) 



V i,s ^ q 



where we have introduced Nambu row- and column- spinors 



ct - fJ 



10\ /01\ /0-i\ /lO 



and 



= 1 1 ; ' ^1 = U j ' ""^ = I i j ' ^3 = V -1 

are usual Pauli matrices. 



(2.11) 
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To derive equations of strong coupling theory, or so-called Eliashberg-type equa- 
tions,c3 we use the following trick. We construct a functional iy[G', D] of Luttinger- 
Ward type.cZI As is known, the functional derivative of such a functional with respect 
to one-particle Green's function G give the self-energy Y,[G,D] (or mass operator), 
as functionals of G and D, and, along with Dyson's equation for Green's function 
G, this procedure leads to a set of self-consistent equations on 12[G,D] and Gl5 
Within the Frolich-type model being used here, "bare" phonon Green's function D 



is fully determined by H^-ph Eq- (2_^)j a-^d coinsidcs with the renormalized one. 



We suppose also that the functional T4^'^[G] determining certain approximate self- 



energy £^[6] for the alloy part of the problem, Eqs. (gj,gjj2j) (i.e., without EPI), 
is known and what remains is to find the contribution of EPI to 14^ [G, D] up to the 
lowest order in EPI coupling with account of those renormalizations of EPI vertices 
by disorder which are consistent with given approximate I]'^[G] (or, which is the 
same, with the approximation to iy'^[G]). 

More formally, we seek for the set of equations of the following general form 

w[g,d]^w''[g] + w''^p''[g,d], (2.12) 
j:[g,d] = j:^[G] + j:^-p^'[g,d], (2.13) 

(2.15) 

G = Go + GoS[G,£']G, (2.16) 

where WIG] and W''^p''[G,D] arc the contributions to the functional T4^[G, D] 
which are due to disorder scattering and EPI respectively; E"^[G] and 'S'^^'p^[G, D], 
the corresponding contributions to the self-energy S[G, D], and in ( 2.16| ) Go is 



"bare" one-particle Green's function determined by non-random part of Hg Eq. ( p.2| 
To construct the functional W'^^p'^ [G, D] we note that to lowest order in EPI 



SD^ 

SW-P'^IG, D] 



5DT 

where n[G, D] is the polarization operator, and 

(5n[G, D 



(2.17) 



5DT 



0, (2.18) 



so that n[G,D] is a functional of G only, H = n[G]. The functional PF'=-p''[G, £>] 
may then be writen as 

W'''P^[G,D]^]^^]yDn[G] (2.19) 

in some formal notation. 
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Now we may give the definition of the renormahzations of EPI which are consis- 
tent with the approximtion for I]'^[G']. We will call the contributions to Ti^~p^\G, D] 
consistent with the approximation to S'^[G] if this contributions are derived from 
the functional ( 2.19| ) where the polarization operator n[G] is consistent with the 
approximation for E*^ [G] in Boym-Kadanoff sence£j'E3 

To obtain explicite expression we devide W^~^'^ [G, D] into two parts 



W 



e—ph 



[D,G] = W"[D,G]+W^[D,G], 



(2.20) 



where [G, D] and [G, D] are Hartree- and Fock-type contributions to the 
generating functional respectively. 

Consider Fock-type contribution first. Contributions to the polarization opera- 
tor 

SW^[D,G] 



W [G] = 2- 



(2.21) 



are shown in Fig. 1. Restoreing the generating functional W-^[G,D] with the use 





r 



Fig. 1. Fock-type contributions to the polarization operator. Black circle denotes "bare" electron- 
phonon three-leg vertex; bold line, the renormalized electron Green's function; the square labeled 
with r, full disorder four-leg vertex. Arrows points the direction of index arrangement. 



of Eq.(^l9|) yields 

W^[G,D] = ~-^ J2 D^y{ips-ipr,^)M2a'{Ga0{Wn^)Gp>a'{Ws) + 

ips ,ipm 

GatJ,{Wm)G^'a'{Ws)Tfj.,yyfj.'{ipm,'iPs)Gp>^'{ips)G^p{ip )}mJ,; (2.22) 

and we depicted the functional W^[G,D] in Fig. 2. Here and in what follows 





Fig. 2. Fock-type contributions to the generating functional. Bold wavy line denotes "bare" 
phonon Green's function and meaning of the other graphical elements is the same as in Fig. 1. 
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implicit summation over doubly repeated Greek indices is implied, 

Ka' = [^3].„.^, (2.23) 

is bare EPI vertex, and we have introduced multy-indices encapsulating site and 



spinor indices (e.g. a, (3, etc. in (2.22)) or site and branch indices (7, 7' in (2.22)). 

The quantity T ajj-^p' a' iyPm^ '^Ps) is full disorder four- leg vertex which obeys the 
Bethe-Salpeter equation 

^a0■0'a'i^Ps,''^Pm) = Uap-p' a' {SPs,'^Pm) + 
+ Uafj.-p.'a'iiPsAPm)Gfj.i,iiPs)Gu'tj.'iiprn)Tul3:l3'u'iiPs,iPrn), (2.24) 

and U^^.y ^'{ips,\pm) is irreducible disorder four-leg vertex generated by the ap- 
proximation used for S*^ [G\ (in fact this vertex is the second variational derivative 
of the functional [G] with respect to G) . 

Apart from the standard Fock-type contribution, there is still a contribution 
of HartreeJype. This contribution is usually neglected following MuttaliJi^iid 
Anderson.lla Howevere as was demonstrated in a series of papers by BelitztSOlla 
this contribution is important and do play role in final expression for Tc. 

To take Hartree-type contribution into account we introduce, in an analogy with 
Fock case, a polarization operator 11^ [G] as 



^ 5W"[D,G] 

(see Fig. 3 for corresponding contributions). We then have 



n^[G]^2 "'^^^^^'"^ (2.25) 



O O <IiI> 

Fig. 3. Hartree-type contributions to the polarization operator. 

W[G,D]^^ D^Y{Q)Ml^.[G^.^{ip^)Gpp.{ip,) + 

{ip„i,ips)Gi3f,' {\ps)Gu'i}' {ips )}a/J;, (2.26) 

and the functional W^[G, D] is shown in Fig. 4. 

Putting everything together one obtains for W'^^'p^[G,D] the expression 

WP'^[G,D] = ^ J2 i?77'(0)Ml'{Ga'4iP™)G/3/3'(ip.) + 

GQ'p(ip,„)G^Q(ip„i)r^t.;i.'p'(ipTO,ips)G/3^/(ips)Gi.'/3'(ips)|Mj,^ 

--^ ^ D^Yi^Ps -^Prn)M2a,Ga'l3'{W7n)Gi3a{iPs)M^,l^{ipm,iPs) (2.27) 
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Fig. 4. Hartree-type contributions to the generating functional. 

and we have introduced renormalized by the disorder EPI vertex: 

Ml,(ip.,ip™) = AC + 

M2^,{ips,ipm) = M2^,+ 
Ual3;l3'a'iiPs,iPm.)Gl3f_c{ips)Gfj,'f3'{ipm)M'^^,{ips,ipm)- 

The use of ( [2.12 - 0^ ) then yields the following set of coupled equations: 

Gap(ips) = Glp{ips) + G'l^(ips)^^i,{ips)G^p{ips), 
{Go = ^PsToSij - (tijTo ~ fj,S,j + eA5ij)T3, 

E^,(ip,) = S^,(ip,) + Spp''(ip,), 
SWIG] 



5Guti{ips) 



+ 



(2.28) 
(2.29) 

(2.30) 
(2.31) 
(2.32) 

(2.33) 
(2.34) 



+r^M';/5/3'(iPs>iPm)Gp'Q(iPs)M2a'GQ/3(iPm)G/3'^'(ipm)A/J,^ + 

+ A/^]^,,G^,/Q(ip„i)G'Q'/3(iPm)AfJ^/G/3'i/'(ips)raa';!.'i/(iPm,iPs) + 
+ 9r^^/ ;;3/3' ; 1/ V (iPs , iPm , iPs ) X 
xGi,>aiiPs)M2^,Gc,>piipm)Gf3'siip,n)M]g,Gs'^'iips)], (2.35) 

E^f,(ip,) = TY,D^Ym{Mj^G^0{ipm)M;l + 

+r^M';/3'/3(iPs,iPm)G/3a(ipm)M2„,Ga'/3'(iPm)G^'i,'(ips)Mj,^ + 
+Mj^,Gpv(ips)Ga/3(ip™)M2„,G/3'a'(iPm)r^'^;/5;3'(iPs,iPm) + 

+c?r^/i';/3/3';'>v(iPs, iPm, iPs)x 
xG^-„(ip,)M2„,Gc«v'(iPs)G/3.5(ip™)Af25',G5.;3(ip„)}, (2.36) 
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where dT denotes the functional derivative of four-leg vertex 



(2.37) 



Note, that in this expression and in similar expressions in what follows the derivative 
is rather a partial variational derivative. A possibility to use partial variational 
derivatives stems from the fact that both fermionic loops and "crossing" propagation 
lines are absent in disorder vertices for the case of quenched disorder. Therefore, 
in the right hand side of the expression ( ^!37|) the arguments ips and ipm label 
two independent functional variables and should be formally considered as being 
distinct. Latter on we will exploit this fact without mentioning explicitly. 



The use of the Bethe-Salpeter equation ( 2.24 ) and definition ( 2.37 ) gives the 
equation on dV: 

dT^^>.ai3.^^>^{ips,ipmAPs) = U^^'.ai3.y;,{ips,ipm,iPs) + 
+ U^^';a0'-^„iips,T^Pm, T^Ps)G jS' s{iPm)G^' ^{ips)T S/S-m' d^Pm, T^Ps) + 



and 



+ C/^^';Q5(iPs,iPm)G55'(iPm)r5'/3;y'^(ip 
-Us' fj.' -aa' iWs,'^P7n)Gs5' i^Ps)Ga' P' {'^Pm)dT fj_S-p^p.^^,^{ips,ipniAPs), 

Uf^f^'^ap-.u'i^iiPs^Pm^Ps) = SU^' ^'..affiipsAPm) / 5G^^{ips) 



(2.38) 



(2.39) 



is irreducible six- vertex due to disorder scattering. Using Eq. ( 2.24 ) once again we 
may express dT through F, U4 and Uq only and after straightforward but lengthy 
manipulations we arrive at the result 

9r^^';a/3;i/v(iPs,iPm,iPs) = U^^' -ap-u' i^iWs , Wm, Ws) + 
+T5'f,':aa' (iPs , ipm)G5S' (jPs)Ga'f3' (iPm ) Uf_,S;f3' P;,y' ^{^Ps , iPm APs) + 

+ t^MM';a/3';Ciy(iPs, iPm, Ws)G 0' s{^Pm)G^' dips)T spy ^ipm, Ws) + 
+r5'^';aa'(iP5-, iPm)GM' (iPs)GQ'/3'(ipm)C^M5;/3'C;»j',y(iP6-, iPm, iPs)x 

xG^^' {ipm)Grfn' {^Ps)T^' p-.yni^P 



■.as{ips,ipm)Gss'iipm)Ts'p-yu{ipm,ips)- 



(2.40) 



Then substitution of ( |2.40| ) into ( p3^j2.36| ) with the use of ( p^j2.29| ) leads to final 
expression for Y^'^^p''-[G, D] (see Figs. 5 and 6) 



spp''(ip,) = s^,(ip,) + s;,(ip,). 



(2.41) 



S^i^fe) = -T^D^yiips 



1P« 



+ ^^MM';<5C;!.v(iPs,iPm,iPs)Gp'a(ips)M2„,(ips,iPm)Ga'5(ipTO)x 
xG^f3{ipm)Mjij,{ipm,ips)Gp'^'{ips)^ 



(2.42) 
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S^, (ip, ) = T ^ D^y (0) { A/;, G„;3 (ip™ + 

+rm«A'^5'5(iPs'iPm)G5/3(iPm)Mj0'G/3'5(iPm)G^'i/'(iPs)Mj,^ + 

(ips,ipm)x 

XGf_t>^{ips)M^^,Gs>a'{^Pni)Ge5i^Prn)Gae'{'^PTn)M^,fi^'t,>{ips) + 
+ Uf^^'.al3;i^'iyi^Ps,iPniAPs)'X 

'xGfj,'s{Ws)M]g,Gs'u'{iPs)Gi3i3'{ipm)Mj,^,Ga'a{iPm) + 
+ ^^AiC;e/3;i.v(iPs,iPm,iPs)Qc'(iPs)G<:'e(ipm)r^'^';Q£'(ips,iPm)x 

xGp'5(ips)M75,G5'iy'(ips)G/3/3'(ipm)M^/Q,/GQ'Q(ipm) + 
+ C^A'M';ae;Cy(iPs,iPm,iPs)G^'^(ips)Gcc'(ip,„)r,,/^/;c'0(ip5,ip„)x 

xG^/5(ip5)M75/G5/,,/(ips)G/3;3/(ip„)M^,^,GQ/Q(ip„) + 

+ C^AtC;«'C':P'^(iPs,iPm,iPs)G^^'(ips)Gee'(ip„)r^'p':Q£(ips,iPm)x 
xr5>';ae'(iPs,iPm)Gp'p(ips)G(;'f (ip™)x 

X Gp' 5 {ips ) M^j, Ga- {ips ) G^j/j' (ip™ G^- „ (ip™ ) + } . (2 .43) 



The formulas (|3o[-|3|) and ( |2ll| - pl3| ) are those which constitute closed set 





Fig. 5. Fock-type contributions to the electron-phonon self-energy. Black triangle denotes the 
renormalized electron-phonon three-leg vertex; dashed hexagon, the irreducible disorder six-leg 
vertex. The other graphical elements are the same as in previous Figures. 



of self-consistent equations in the case of EPI renormalizations taken into account 
consistently with the approximate S^[G] given. 

To obtain the equations for the superconductive transition temperature Tc one 
should linearize the equation for S[G, D] in anomalous part which is assumed to be 
small near Tc. This is formally achieved by expanding S[G] in a functional series 
in SG, 5G being anomalous part of G, and then using Dyson's equation to express 
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Fig. 6. Hartree-type contributions to the electron-phonon self-energy. 
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SG back through 6T,, where 5T, is the anomalous contribution to S[G, D]. We have 



6j:[g,d] 



SG 



6G = 



nGn,D] 



6G 



G=G„ 



(2.44) 



G=G„ 



(subscript "n" means that the corresponding quantity is taken in the normal state of 
superconductor) and equating contributions of order 0(5E) in left- and right-hand 



sides of Eq. ( 2.44 ) we obtain formal equation 

(5S[G, D] 



SG 



GnS'YGn, 



(2.45) 



G=G„ 



and Tc may be found from solubility condition to this equations. Note, that the 
solubility condition may be written as equality of some formal determinant to zero, 
but we omit explicit formal expression of such a form. The diagrammatic represen- 
tation of ( ^.45 ) in the case under discussion contains ninty seven terms and is too 
cumbersome to be written down here. 



Hence, normal part of Y,[G,D] near is still determined by Eqs. (2.32 2.33) 
and ( pll}fl3| ) with one-particle Green function G taken in normal phase of su- 
perconductor, and we have formal equation ( p. 45 ), from which the expression for 
Tc may be derived after suitably fixing approximation for E'^[G]. 

In the conclusion of this Section some remarks are to the point. 

Obviously, the approximation to I][G, Z?] given by ( 2.32 , ^.33 , 2. 41 - 2. 43| ) is con- 
serving in a whole by construction. 

The expression for [G, D] ( 2.27 ) does not actually depend on particular 

alloy model since the assumptions on the alloy model enter into this expression 
only implicitly through irreducible four-leg vertex, so that the form of W'^~p^[G, D] 
is independent of the disorder type. We will demonstrate the correctness of the 
statement above in subsequent Section for the substitutional disorder of a general 
type carrying out configurational averaging in a formally exact way. 



The equations Eqs. ( 2.32 , 2.33 ) and ( 2.41 - 2.43 ) enable also to consider "incon- 
sistent" renormalizations of EPI vertices. In this case we assume that E*^[G] is 
calculated within certain conserving approximation and renormalizations of EPI 
are taken into account within another, but still conserving approximation for S'^[G]. 
This happens, for instance, when discussing the so-called Anderson's theorem where 
one assumes that disorder does not influence the EPI at all. The last, in turn, im- 
plies that the approximation for Y.'^ [G] within which we renormalize EPI vertices is 
E'^[G] = 0, and therefore irreducible four- leg disorder vertex and of course higher 
vertices vanishes identically. 

3. T Matrix Approach to the Superconductivity of Disordered Alloys 

In this Section we derive the equations ( ^.30 - 2^3^ ) and ( 2.41 - 2^4^ ) for the case of for- 



mally exact configurational (disorder) averaging. Making use of multiple scattering 
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theory, or so called T-matrix approach, we establish relations with the approach 
of the previous Section and demonstrate that the Eqs. ( ^.3C - 2T3^ ) and (2.41-2.43) 



are actually the case for general type of substitutional disorder when the disorder 
averaging is carried out exactly. However, if the matter is a specific approximation 
within certain model of substitutional alloy, it is more convenient to formulate the 
approximation in terms of generating functional, so that the approaches of Section 2 
and of this Section complement each other. 

For the first time the application of T-matrix approach to the problem of alloy 
superconductivity has been given in Refs.i'§0 where the model of a binary alloy 
was considered with local (single-site) BCS-type interaction and, hence, with purely 
local superconductive pairing. In these papers the equation for the self-energy was 
derived within the framework of Coherent Potential Approximation JCPA) and 
resulting BCS coupling renormalized by disorder was presented in Ref.Ll. 

Theresulta^of this Section are a strightforward generalization of the approach 
of Refsfl'EHaBlia to the case of non-local superconductive pairing owing to EPI with 
bare EPI vertices being independent of disorder. 

For particular alloy configuration, consider one-particle Matsubara Green's func- 
tion g(T) with matrix elements 

{gML/3 = -<T.c„WcJ(o)>ff' (3-1) 



where < ... >h means usual Gibbs averaging, with H being the Hamiltonian ( |2.l| ). 

After transforming to Matsubara frequencies and taking EPI into account one 
has for g~^(ips) the expression 

g-i(ip,) ^ ipsTo - Ht3 - V'^-P'^iips), (3.2) 



with H being the matrix with elements (2.3); V*^ ^''(ips) being the effective inter 



action due to EPI in this disorder configuration, possessing functional dependence 
on g(ips)- 

To derive equations required it is sufficient to know the functional V'^^P''(ips) 
to the lowest order in EPI of the self-consistent perturbation theory. To this order 
the expression for V^^P^{ips) has the form 

V--P''{ips) = V"{ips) + V^iips), (3.3) 



V^piips) = TY,D-n'im2f}go.'f3'iiPm)Mj:^,, (3.4) 

Vf/jliP^) = -T^D^yiips - \Pm)Ml^,ga'f3'(\Pm)Mj,p, (3.5) 
iPm 

where V^(ips) and 'V^{\ps) are Hartree and Fock contributions to the effective 
potential respectively. 
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For the con figuration a] avR raging to be carried out we introduce T-matrix. In 
the usual manner one hasE2rEZI 

T(ip,) - (V(ip,) - S(ip,)) + {V{ips) - S(ip,)) G{ips)T{ips), (3.6) 
T(ip.) = (V(ip,) - S(ip,)) + T(ip,)G(ip,) (V(ip,) - S(ip,)) , (3.7) 



where 



V(ip,)=V^+V^-P^(ip,), 



(3.8) 



V*^ is random contribution to Hamiltonian ( |2.2| ), depending on disorder configura- 
tion, and 

Siips) = G(ip,) + G(ip,)T(ip,)G(ip,), (3.9) 
< g(ip,) > = G(ip,) + G(ip,) < T(ip,) > G(ip,), (3.10) 

G(ip3) - Go(ip.) + Go(ip,)S(ip,)G(ip,). (3.11) 

Here < . . . > stands for configurational averaging with particular probability dis- 
tribution and Go(ips) is completely determined by non-random part of H ( ^.2|) . 
For self-consistency S(ips) is fixed by the condition 



< T(ip,) 0, 
and, hence, from ( ^^9 3.1(]| ) and ( |3.12| ) it follows 

< g{ips) >= G{ips). 



(3.12) 



(3.13) 



To carry out configurational averaging we must use two apparently different 
schemes of averaging: one for the alloy part of the problem and another for the 
effective scattering owing to EPI (the latter scheme is reminiscent of Virtual Crystal 
Approximation (VGA) for ordinary "alloy" problems). To this end we decompose 
S(ips) according to natural partitioning of V(ips) into two parts owing to disorder 
scattering and EPI respectively 



S(ip,) - S^(ip,) + s^-p"(ip,). 



and then rewrite Eqs. (3.6.3.7) for T-matrix in the following form 

T(ip,) = Ti{ips) + T2{ips), 
Ti(ip3) = Vi(ip,)+Vi(ip,)G(ip,)T(ip,), 
T2(ip.) = V2(ip,)+V2(ip.)G(ip,)T(ip,), 



where 



Vi(ip,) = Y'^-'E^ip,), 

V2{ip.) - V'=-P''(ip,)-S'^-P"(ip,). 



(3.14) 



(3.15) 
(3.16) 
(3.17) 



(3.18) 
(3.19) 
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After introducing partial T-matrices T*^ and T*^ by 



T^fe) - Vi(ip,) + Vi(ip,)Gfe)T^(ip,), 
T^-P''(ip,) - V2(ip.) + V2(ip.)G(ip.)T^-P''(ip,), 



equations (3.16,3.17) for Ti and T2 transform as 

Tife) = T^fe) + T^(ip,)G(ip,)T2fe), 
T2(ip.) = T^-f''(ip.) + T^-P'^(ip,)G(ip,)Ti(ip,). 



(3.20) 
(3.21) 



(3.22) 
(3.23) 



Bearing in mind the lowest order in EPI let us find the full T-matrix (3.6, p77]3.15| ) 
up to the first order in V2. In this case Eq. ( |3.2l| ) becomes 



T^-P^tip,) =V2(ip,). 

To this order in V2 the result for the full T-niatrix reads 

Tixps) = T'iips) + V2(ip,) + V2(ip,)G(ip,)T^(ip,) 
+T^(ip,)G(ip,)V2(ip.) + T^(ip,)G(ip,)V2(ip,)G(ip,)T'=(ip,). 

Fixing up to now arbitrary quantity S'^(ips) by the condition 

<T'=(ip,) >=0, 



(3.24) 



(3.25) 



(3.26) 



we obtain from ( 3.25 ) the equation for P^iips): 

S'=-P''(ips) =< Y^-P^{\ps) > + 
< V^-P''(ip,)G(ip,)T^(ip,) > + < T'=(ip,)G(ip,)V^-P''(ip,) > 
+ < T^(ip,)G(ip,)V2(ip,)G(ip,)T^(ip,) >, 



where the use was made of Eqs. (3.1 2B^9| ) and ( [3.26| ) 
Introducing in a standard fashior 



10 full four- and six-leg vertices 
rpM';i"''(iPsjip™) = < r^^/(ips)T;^'(iPm) >, 

^aa';l3p'-n7'{Ws,iP7n,iPs) = < T^^, {ips)Tp^, {ip„i)T^y (iPs) > , 



(3.27) 



(3.28) 
(3.29) 



writing down indices (the summation over doubly repeated Greek indices is implic- 
itly meant here and in what follows) and then substituting (3^,3^) and ( ^.S| , 3.10| ) 
into ( 3.27 ) we obtain 



S^,.''''(iPs) +r^M';!^'>'(iPs)(iPs)G'A"a(iPs)S^/ {ips)Gp^{ips) = 



4i^(ip,) + 4i^(ip,), 



(3.30) 



where 



^i'K^Ps) = -TY,D,y{ips " ip,n){M;^Go.0{ip^)Mjl+ 
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and 



M 



;/3/3'(iPs,iPm)Gp'Q(ips)*^2a'^a'/3(iPm)G;5',,'(ip„)Mj,^ + 

(ips,iPm) + 
fe)} 



E(2(ip,) = T^i?^y(0){Mj„G„MW)M;:- 



(3.31) 



),'u'{iPs)Ga[3{ipm)M2a,Gi3'a'{iPm)^u'u;f3P'{iPs,'^Prn)' 
xGp'a(ips)Af2a'GQ'^'(ips)G/3'5(ipm)M2'5,G5'/3(ipm)|. 



(3.32) 



The equations (3.30-3.32) are essentially a generalization of corresponding equa- 
tions of RefsB00 to the case with EPI-induced superconductive pairing and when 
configurational averaging is carried out exactly, and they may be used to discuss 
the phenomenon of the superconductivity of alloys on the grounds of certain ap- 



proximate solution for ( 3.20 ) and S'^(ips) ( 3.26 ). In particular, for the case 
of diagonal alloy disorder the replacing T*^ with t in ( 3.3G| - 3.32 ), where t is the 
single-site T-matrix, 



t{ips) = {e^T3 - T.'^iips)) + (£jT3 

3U.ti.ons 
wit 



i:'{ips))Gu{iPs)t{ips), 



(3.33) 



retaining single-site EPI contributions only and using CPA expression for E*^, one 
may obtain the equations of Refs.Btl with the exception that superconductive pairing 
is due to EPI, not BCS interaction. 

However we prefer the equations for the self-energy S(ips) in the form, which 
does not rely upon T-matrix and gives a possibility to analyze arbitrary conserving 
single-site approximation for S(ips) within an unified scheme. Therefore we rederive 
the expressions ( |2.4lH2.43| ) for i:^-P^{ips) from ( |3.30H3.32| ) to demonstrate that the 
two expressions are in fact equivalent in the case of exact configurational averaging. 

In a standard way, the irxeducible disorder four-leg vertex, U4, is introduced by 
the Bethe-Salpeter equationcSlO 



Tf^^'-„'^{ips,iPm) = Ui_,^>.y^{ips,ipm) 
{ips, iPni)Gaa' {w s)G piipra)^ a' fj^' mu' P' {Ws ^Pn 



(3.34) 



(note the inverse order of introducing of full and irreducible vertices in this and 
in the previous Sections). Multiplying ( ^.30 ) by the combination GU4G (where 
indices must be identified explicitly) we can express < T'^GS'^~^'''GT'^ > in the 
left-hand side of ( ^.3C| ) through U4, Fg and G. Substituting the result obtained 
back into ( 3.30|) and collecting terms we then arrive to the result which reads 



SP'"'(ip.) = E^,(ip,) + S^,(ip,), 



(3.35) 
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Mll^,Gfj^'aiiPm)Ga'piipm)Mjfj,Gi3>^'{ips) 

;/3/3' ; 1^ V (iPs , iPm , iPs ) X 
xG^,a.iiPs)M2^,Ga,'p{ipm)G^,siiPm)M]^,Gs','iips)} (3.36) 

S^.fe) = Tj2D^YiO){Mj^GMiPm)M;ii + 

Tf,f,'.f3'f3iips,iPm)G0a{iPm)M2^,Ga'/3'{ipm)G^'ty'{ips)M2,^ + 
M2,,'G,,'^'{iPs)Ga0{ip,n)M2a,G(3'a'{iPrn)'^i^'i^-pp'{iPsAPrn) + 

;/3/3' (ip* ' : ) x 
xG^.a(ip,)Af2„.G„'.'(ip.)G;5.5(ip„)M75',G5.;5(ip™)}, (3.37) 



and 



Qtiti';Si;'^'ui^P.s,iPm,^P.s) = r^M';5«:;'';'(iPs'iPni'iPs)~ 

Uf_,^-^':,{ips:iPs)G^^'{ips)T^^^>.M'iipsAPrn)~ 

m J iPs j 

'^MC;«''(iPs'iPs)G'«'(iPs)G'KK'(iPs)r4V';55';i/K(iPs,iPm,iPs). (3.38) 

Here the expression (3.38) is nothing but a variational derivative dT of r4 ( 3.28| ) 
with respect to G. Indeed, the variation of the T-matrix ( 3.20| ) with respect to G 
gives 

= -T" |(5 ((V^ - S'^)"^ - /(5g| T'^ 

= T-^ |- (V" - S-^)^^ (5SV5G [Y^ - S-^)"^ + SG/SG^ 
= - (1 + T-^G) (5SV(5G (1 + GT-^) + T'^(5G/(5GT% 

and 



SGa^ 



SGjs 



(3.39) 
(3.40) 



Introducing 



C^Q7;<5/3(iPs,iPs), 



(3.41) 



SGjsiips) 

expHcitly writing down indices in ( 3.39| ) and then substituting (3.39) into the fol- 
lowing expression 



l9r^;^';55';i/'i/(iPs,iPm,iPs) =< 



ST^-^^-jips) 
SG^^{ips) 



TsS'i^Pm) >, 



(3.42) 
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one obtains 



dr 



U^,!=-yt,{ips,ips)G^^'{ips)r^>^,.ss'{ipsAp,n)'- 
Uf^f,' ■^',y{ips,iPs)G^^' {ips)Tss' -u'^i^p 



where we have used ( p.4lD and definitions ( |3.28| ) and ( |3.29D 



(3.43) 

The four-leg vertex 

(3.41) is completely determined by the Bethe-Salpeter equation and, as shown in 
Appendix A, the definition (3.41) is in fact the Ward's identity. 



Comparing Eqs. (3.43) and (3.3J) one sees that they coincide with each other 



and, consequently, (3.35-3.37) and (2.34- pT36 ) are identical. In Appendix A we will 



demonstrate that Eg. ( 3.26 ) determining the self-energy 'E^{ips) can be represented 



in the form (2.33). Simple repetition of the transformations of the previous Section 



(Eqs. (^.34| - |2.43[) ) then leads to Eqs. (|2.4l|-|2^) . This completes the demonstration 



(2.33 


,^.41-2.4S) and ( 


3.26, 


3.30 


-3.32) 



the case of exact configurational averaging. 

Note that the case of exact configurational averaging we dealt with in this Section 
may be considered as an extreme case of conserving approximation in the sense 
that S'^[G] coincides with the exact value, and the condition of (/)-derivability of the 
approximatioriEj 

= (3 44) 

SGjsiips) 6Gi3a{ips)' 

can be readly verified (see Appendix A). Therefore, as for conserving approxima- 
tions, formal results which will be obtained in the rest of the paper are valid for 
this case too, and we will not discuss it separately in the following. 

Concluding this SectiQrLjfit^us stress some important points. 

In contrast to Refs.Hl3fl&t3, where site-diagonal contributions to the effective 
scattering potential (which come from BCS-like interaction) have been treated on 
an equal footing with disorder contributions, that is within the same approximation 
(in fact CPA) , we rather apply the same Virtual Crystal type averaging procedure 
to both site-diagonal and site-off-diagonal EPI contributions to V(ips)- The decom- 
position of T,[G, D], Eq. (|t|), and, hence, of the fuU T-matrix, Eqs. (|T5|-|t|), 
gives a possibility tocarry out two different averaging procedures. Thejjery con- 
tent of the Migdal-Eliashberg theory, based on the addiabaticity theorem,C3 dictates 
such the way of partitioning the self-energy D] into the parts owing to disorder 
scattering and EPI, as well as it dictates the virtual crystal type averaging proce- 
dure for the EPI contributions. Then the Eliashberg-type equations acquire the 
character of effective field equations with two effective fields I]'^[G'] and T,'^~p^[G, D] 
coupled via Dyson's equation for the one-particle electron Green's function G. 

Suppose for a moment that we have decomposed the full self-energy into, for 
example, three parts owing to disorder scattering, site-diagonal and site-off-diagonal 
EPI contributions, respectively. Then, if we use CPA-like, or even exact, averaging 
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scheme for site-diagonal EPI contributions, this will give results which are formally 
beyond the accuracy of the Migdal-Eliashberg theory as they contain, besides the 
lowest order contribution, higher order EPI contributions to Y,'^^p^[G,D]. 

Another important feature is that the decomposition (|3.8|) of V(ips) and S[G, D], 



Eq. (3.14), is the only possible for both the Ward's identity (p. 41) and (/)-derivability 



criterion (3.44) for E'^[G] to hold. Rather different approaches of this Section and 
of Section 2 are then in tight analogy with each other: the form of the equations are 
the same and detailed structure of underlying quantities may be established using 
the connections between disorder vertices and the T-matrix. 

4. Ward's Identities and the Anderson's Theorem 

As a first application of the approach being developed in this paper we consider the 
so-called Anderson's theorem which concerns the influence of nonmagnetic disorder 
on the superconductive transition temperature T^- 

This theorem0@0 was established within BCS model and Virtual Crystal Ap- 
proximation (VGA) for dilute nonmagnetic impurities^ and for concentrated non- 
magnetin impurities taken into account within Coherent Potential Approximation 
(CPA).E3 For the latter case the theorem was discussed on more general grounds.Q 

For isotropic superconductor where pairing owing to EPI is primarily of s- 
character there is a standard consideration in the literature 3 of how dilute non- 
magnetic impurities (essentially the case of weak alloy disordering) influence Tc. 

For the case where disorder scattering strength is arbitrary, which was treated 
by using Average T-matrix Approximation (the concentration of the impurities is 
small none the less), and in the absence of EPI renormalizations this question was 
discussed in connection with one-particle state density variations near the Fermi 
energyO It was demonstrated that, while usual cancellation of anomalous disorder- 
scattering contributions holds in the equation for the gap function, the asymmetry 
of the normal-phase one-particle state density near the chemical potential leads 
to non-zero (in fact, non-constant) spectral shifts. In this case, equations for the 
gap function and for the spectral weight renormalization function Z{ips) no longer 
decouple which results in the disorder contributions to the gap function via explicit 
dependence of the corresponding equation on Z{ips). Then, as was concluded,EJ in 
the presence of any structure of the density of states (DOS) near the Fermi level, the 
Anderson's theorem breaks down, which seems to be in contrast with the situation 
for the BCS model. This point needs further clariflcation. 

Further we shall see that the account of the Fermi-surface energy dependence 
and its electron-hole symmetry is not directly related to the Anderson's theorem. 
It is strong k-dependence of the superconducting order parameter (Aj^(ips) is not 
constant over the Fermi surface) only which may lead to the suppression of super- 
conductivity by non-magnetic impurities in the case of weak disorder. However from 
the poit of view of the standard considerationsQ this situation does not correspond 
to s-pairing (basis function in the expansion in partial waves for an anisotropic 
Fermi surface is unity for s-channel!). So, the term "s-pairing" will be understand 
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hereafter as follows: the gap function Aj^ has the same point symmetry as the 
"bare" electronic spectrum ej^, that is Aj^ = A(ej^). One may expect that this holds 
true for all conventional superconductors. 

Within the Bogolyubov's formulation of the BCS-likc model, that is for the 
model with disorder-independent attractive four-fermionic interaction which is al- 
most of infinite range and has a separable kernel, and after adding disorder-de- 
pendent two-fermionic contributions, the Anderson's theorem may be understand 
basing on qualitative considerations. It states that the expression for Tc has usual 
BCS form with the normal state one-particle spectral density renormalized by the 
disordem replacing bare one. Indeed, by analogy with the case where diorder is 
absent ,e3 one can construct an approximate Hamiltonian of effective- field type, 
which gives asymptotically exact (in the thermodynamical limit) solutions for the 
thermodynamical quantities and for correlation functions as well. For fixed disorder 
configuration the effective field (in fact, the superconductive gap) which enters the 
approximate Hamiltonian, is spatially homogeneous and can be determined through 
the self-consistency conditions containing anomalous averages of fermionic operators 
forming Cooper pairs in this disorder configuration. The anomalous averages are 
then expressed in terms of the one-particle spectral density which is non-negative 
and is the same for both spin-up and spin-down electron states because, as was 
pointed out by Anderson, the states of electrons forming Cooper pair are time- 
reversal images of each other. Hence, the effective field may depend on disorder 
via the corresponding dependence of the spectral density, but the "disorder" fluc- 
tuations of the effective field are strongly suppressed in the thermodynamical limit 
because of infinite range of the interaction, so that the effective field can be safely 
replaced by its configurationally averaged value in the self-consistency equation. 
Residual disorder averaging touches only the spectral density; the self-consistency 
equations become similar to those of usual BCS theory; however they contain the 
renormalized spectral density instead of bare one. The statement of the Anderson's 
theorem then follows. Note, that for this particular model the conclusion takes 
place irrespectively of the strength and type of quenched substitutional disorder 
and also of whether renormalized and/or bare state densities possess any structure. 

The situation for the model (2T) is somewhat different. In this case effective 
EPI-induced electron-electron interaction is rather of short or intermediate range 
than of (infinitely) long range, "disorder" fluctuations of the electron-phonon part 
of the self-energy are not small in general and lead to renormalizations of "bare" 
EPI vertices. But if one neglects completely these renormalizations the statement 
may be established which is quite analogous to the Anderson's theorem for the 
BCS-fike model. 

In this Section we demonstrate that, for superconductive alloy with isotropic 
s-type pairing owing to EPI and for arbitrary conserving approximation for the 
disorder self-energy S'^[G], quenched substitutional disorder influences supercon- 
ductive transition temperature Tc only through normal-state one-particle spectral 
density renormalized by all the interactions in the system, provided disorder does 
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not change "bare" EPI vertices. Surely, this makes the theorem to be rather a for- 
mal statement because the effects of EPI- vertex renormalizations, as we will see, are 
not small except for the case of weak alloy disordering (dilute alloys). Note, also, the 
Anderson's theorem is not a statement of mathematical regour but holds to certain 
approximations usually accepted in conventional theory of superconductivity, 
^^he most general form for spatial Fourier transforms of S'^(ips) and Y,'^~P^{\ps) 

S|.(ips) = iPs[l-7k(iPs)]'ro + Xk(iPs)T3+(/)|,(ips)ri, (4.1) 
Y^'^^^'^CiPs) - ip.7k(iP3)[l - %(iP3)]^o + Xk"'"(iP^)^3 + C^'(iPs)ri, (4.2) 

where only ti contribution to the anomalous part is retained (in the absence of 
external magnetic fields we can always make T2-contributions vanish by suitable 
gauge fixing because usual gauge invariance under phase transformations is pre- 
served within 0-derivable approximations) and 

7kfe)-l-^^^ ^ ■ (4.3) 

The representation (O) for Y.'^^^iips) differs from the standard onc@ in that the 
disorder contribution have been picked out explicitly into the factor 7j^(ips) (which 
is nothing but the disorder self-energy correction). 




Fig. 7. Standard contribution to the normal electron-phonon self-energy without disorder renor- 
malizations. 




Fig. 8. Full contribution to the anomalous self-energy. Dashed circle represents the anomalous self- 
energy; dashed square, the irreducible disorder four-leg vertex. All the other graphical elements 
are the same as before. 



Linearized Eliashberg-type equations without disorder renormalization of EPI 
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read (see Figs. 7 and 8 for the normal and anomalous contributions to the self- 
energy) 

iPs7k(iPs) [l - ^ki^Ps)] = ~Tc ^ Aijq(ips - ipm)*q(ipm)ips7q(ipm)^q(ipm), 

(4.4) 

X^'''' i^Ps) ^ -Tc ^ Ai^.q(ips - ip™)*q(ip„)£q(ip™), (4.5) 
q,ip™ 

0k(iPs) ^ \qi^Ps - iPm)*q(iPm)'/'q(iPm) + 



■I^t^kq(ip.)*q(ip.)0q(iPs), 



(4.6) 



where 



(4.7) 



is the full anomalous contribution to I]]^(ips), Uy}^{ips) is proper irreducible four-leg 
vertex related to anomalous contribution to I]'^(ips), and 



xV^'^i^Ps)^ 



(4.8) 

ek(iPs) = tk + xli^Ps) + x'^'^iiPs), (4.9) 

with t-^ being bare electron dispersion defined by non-random part of ff*^, Eqs.( |2.2| , 
2.8). Quantities 7^(ips) and x^(ips) describe disorder scattering and satisfy formal 
equation: 

• h r M ^ e,. ^ SW-[G] 



<5Gk(ip.) r^r^ 

Explicit form of this equation depends on the approximation used for ^^[0] or 
S'^[G'] and is not important for further consideration. Note that we have also in- 
cluded Hartree contributions to the renormalization of the chemical potential. 
At last, the expression for Aj^q(ia;„) reads 



f°° 

Akq(iw„) = dQa^F-^^^in)-^-^^ 



(4.11) 



Here a^Fj^q(rJ) coincides with the usual Eliashberg function up to a factor of di- 
mension energy and has the form 



%q(^)-EK<ir{-^-^l;iq 



(n + io) 



Setting in (A. 10) of Appendix A 

G' = G,,-T(ip.), G-G,,a(-ip,), 
S' = S^^.^(ip,), E = S^-,^(-ip,), 



(4.12) 



(4.13) 
(4.14) 
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where all quantities are taken in the normal state, and Fourier transforming one 
then obtains 

S^^fe) - ^'_^yips) = ^kq('^^^ (GqT(iP.) - G_qx(-ip.)) , (4.15) 

q 

or, with account of the definition ( [4.3[ ) for 7j^(ips), 

7k(ip.) = 1 + E C^iq(iPa)^q(ipa)^q(ipa)7qfe), (4.16) 

q 

where the use was made of Eqs. (4.1,1.2.01,4.9) and of the expression 

and it is essential that the irreducible four-leg vertices which enter Eqs. ( |4.6| ) and 



( 4.16 ) are the same. 

Introducing in a standard fashion normal-state one-particle spectral density 

^kaW=--^Ii^Gk,(£ + iO) (4-18) 

and noting that 

^kTW^^-ki(^) (4-19) 

(this is, in fact, a consequence of time reversal symmetry of the problem), we can 
write 

7k(iP.)^kfe)*k(i?'-) = *k(iP«)' (4-20) 

*k(ip.) = / *4t^ (4-21) 
J pi+ e 

and 

Yi,{ips) - e]^{ipsW^iips) = I de'^^. (4.22) 

Here we fix the zero of the energy self-consistently at the true chemical potential. 
The system of the equations (|4.4|-[4.6|) then reads 



^Psly^i^Ps) [1 - Z]^{ips)] = -Tc E Ajj.q(iPs - iPm)*q(iPm)ip.5, (4.23) 

q,ip„ 

"'(iPs) = -Tc ^kq(iP« - iPm)Y]^{^Ps) (4.24) 
q.ip™ 

and 

j]^{ips)Z]^{ips)A]^{ips) ^Tc Y ^kqfe ~ iP™)*q(iP™)^q(iP™)+ 

q,ip™ 

+ E C/j^q(iP.)*q(iP.)^q(iP.)7q(iPs)Aq(ip,), (4.25) 
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where 

Ak(iP.) = ^T^Pl-, (4-26) 

is, as usual, the gap function. 
Solving ( p3| ) for Z]^{ips): 

= 1 + T^7k(iPs)^c ^ •^kq(iP'* ~ ipm)^'q(ipm)ipm (4.27) 



and inserting the result into (4.25) we rewrite the equation for Ajj(ips) in the form 

^(9kq(ips)7q(iPs)Aq(ips) = ^ Ai^q(ip3 - ip™)^'q(ip„)Aq(ip„)- 
q q,ip™ 

_ Ak(iPs) y^ ^ Akq(iPs - ip™)*q(ip™)ip™, (4.28) 

where 

Qkqi^Ps) - '5kq(iP^) - C/iij(ip.)^q(ip.)*q(ip.) (4.29) 

and up to now no approximations has been made. 

Usual way to proceed further is to introduce harmonics which constitute or- 
thogonal and full set at some "bare" Fermi surface, that is so called Fermi surface 
harmonics. We have J0 

Y^Fj{k)Fj,{k)S{e - 4) = No{e)Sjj,, (4.30) 
k 

iVo(e)-E'^(^-4)- (4-31) 
k 

The expansion over Fermi surface harmonics for various k-dependent quantities 
reads 

= E / - el)Fj{k)Cj{e), (4.32) 



^kq = E / dede'Sie- e*^)6ie' - e*q)Fj{k)Fj,{q)Kjj,{ee'), 
J J' 



(4.33) 



^kq = E / - 4)F^(k)i^^,q(e), (4.34) 

7 ^ 



^k^k = E / de5{e-el.)Fj{k)Y,CjLL'AL{e)BL'{e), 

J •' LU 



(4.35) 



E^k^k= / de7Vo(e)E^L(f)^i(e)> (4-36) 

1_ T 
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where coefficients Cjll', being Clebsch-Gordan coefficients for the harmonics, obey 
the property 

CoLL' — Clol' — Cll'q — Sll', (4.37) 

ej^ is "bare" dispersion and we choose e t = tjj., and ij^ is bare electron dispersion 
defined by nonrandom part of H"^, Eq. (2.2). 

With the use of (4.30- O6| ) the equation ( 4.28| ) can be represented as 

f de'No{e')Qjj,{ee':ips)Y,C.j'LL'^L{e';ips)^L'ie';ips) = 

J' •' LL' 

= TcY.Y. J de'No{e')\jj,{e€';ips^ip^)x 
X Cj^LL'^Ue; iprn)AL' (c'; ip^) - T, ^ CjlL' ^^^^^ X 

LL' LL' 

X ^ / de'Noie')\.j'L'i££';ips - ip,„)l^L(e'; iPm)ipm (4.38) 
J' 

Since s-type pairing dominates for the case of isotropic superconductor and, more- 
over, AL(e;ips) is almost independent of e in a narrow energy region |e| < cod the 
following approximation takes place 



AL(e;ips) — > Ao{e*;ips) = A(ips 



(4.39) 



After substituting ([4.39f ) into ([4.38| ) and using ( |4.37| ) we get 



Y <3o,q(e*; ips)7q(iPs)A(iPs) = Tc Y ^o,q(e*; ips - ipm)^'q(ipm)A(ip,„) 

q q,ipn^ 
A(ip, 



iPs 



-Tc Yj •^o,q(e*;iPs - iPm)*q(ipm)ipr, 



q,ip„ 



The Ward's identity ( 4.16[ ) can be rewritten as 

Xl'5kqfe)7q(iPs) = 1, 



(4.40) 



(4.41) 



where Q]i^q{Ws) is defined by ( 4.29|). After carrying out expansion over Fermi 
surface harmonics, the expression ( 4.41 ) becomes 



and thus 



X!Q-'.q(^'iPs)'Tq(iPs) = Sjo, 

q 



X!'9o,q(e,iPs)7q(iPs) = 1- 



(4.42) 



(4.43) 
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The use of (4.43) in the left-hand side of (4.40) immediately yields 



A(ip,) J 1 + V Ao,q(e*: 



ips - ipm)*q(iPm)i?'" 



Tc ^ Ao,q(e*; ips - ipm)*q(ipm.)A(ip„i). 

q.ip„> 



(4.44) 



This final expression clearly shows that disorder contributions enter the gap- 
function equation only through ^'j^(ips) and, thus, through renormalized normal- 
state one-particle spectral density ^^(e) (Eq. ( 4.18 )). Disorder contributions to 
anomalous part of the self-energy cancel the factor ^^^{ips) in the left hand side of 
( 4.23 ) owing to exact Ward's identity ( 4.16 ) and we may conclude that it is a type 
of Ward's cancellations which leads to the Anderson's theorem. 



5. Reduced Equations for the Isotropic superconductor. The Abrikosov's 
Identity 

In this section we will derive a set of linearized Eliashberg-type equations within 
isotropic approximation for the gap function Ajj(ips), that is for the case where 
k-dependence of Ajj(ips) is rather weak and may be completely neglected in the 
narrow energy interval iwu near the true chemical potential. Analogously to the 
case where disorder contributions to EPI matrix elements have been neglected, we 
will demonstrate that, owing to exact Ward's identity, anomalous contributions to 
the self-energy, which come from disorder scattering, may be eliminated from the 
full equation for the gap function. 

With the use of (4.1-0|) linearized Eliashberg-type equations can be formally 
written as 

iPs7k(iPs)[l - ^kfe)] = ^kq(iPs'iP'«)' (5-1) 



-ph 



(ips) = -Tc ^ L 



q,ip 

(x) 
kq 



{ips,ipm), 



(5.2) 



(l^liC^Ps) = Tc ^ L|^^^(ips,iPm)*q(iPm)</!'q(iPm)- 



qap„ 



E ii,'^(ip.,iPm)*q(iPs)0q(iPs) + 5]C/j^ij(iPs)*q(iPs)</'q(iPs), (5.3) 
qip™ q 



where we have denoted contributions to Zj^(ips) and xj^(ips) as Lj^^(ips, ip„j) and 

Tj^q(ips, ipm) in Eqs. (|0|) and ( |5^ ) respectively. Analogously, in ^kq(^?'s' ^^'™) 
means the contributions where anomalous part 4)-^^{ips) is summed over both quasi- 
momentum q and frequency ipm , and L^^{ips,ipm), the contributions which con- 
tain 0jj(ips) summed over the quasimomentum only. Note, bytheway, that while 
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f 2") 

both Hartree- and Fock-type diagrams contribute to ij^^(ips, ipm), only Fock-type 



diagrams give non-zero contribution to _Lj^^(ij3s, ip,„). The last term in (5.3) is 

anomalous disorder scattering contribution to (j)-^(ips), and we have also explicitly 
written down sums over internal quasimomentum and frequency. 



Making use of Eq. ( 4.19 ) we can express ( p.3| ) in terms of the one-particle spectral 
density. A little manipulation yields 

^kq(ip.,iPm)*q(iPs)Aq(ip,)-|-^C/j^^q(ip,)¥q(ip,)Aq(ip,), (5.4) 



q,ip„ 



and 



q 



(5.5) 



Again, the gap function A]^(ips) has been introduced, see Eq. ( 4.26| ). 
Formal solution to the equation (5.1) for Zj^(ips) has the form 



= 1 + -T^^c V L|^^(ips,ip™). 



(5.6) 



The substitution of (|5.6D and (4.16) into (5.4) gives 



[ q,ip™ J 

+Ak(ip,)^C/^yip,)¥q(ip,) = 

q 

= ^ |L^(ips,iPm)*q(ipm)Aq(ip,„) + 
q,ip™ 

+^kq('^^' iPm)^q(ip«)Aq(ips)} + J] C/kq(iPs)*q(ips) Aq(ips) (5.7) 

and, after expanding over Fermi surface harmonics ( 4.3C ) and using 
one obtains 

53 AL(e;ips) < <5lj + CjLL'T^Tc ^ L^''?q(e; ip^, ip„i) > -|- 
ii^' [ q,ip™ J 

+ CjLL'^L{t- IPs) J2 ip,)^q(ip,) = 

LL' q 

= ^-EEE / '^f'^o(e')Cj'LL'X 

ip„, J' LL' 
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X ^^L[pj,{ee'; ips, ipm)*L(e'; ipm)AL/ (e'; ipm)+ 
+L[f], (ee'; ip„ ip„)l^i(e'; ip,) Ai-(e'; ip,)} + 

C/j^, (ee'ip, ) (e' ; ip. ) Ai. (e' ; ip. ) • (5.8) 

7/ T T ' 



J' LU 



As in the previous section, a weak dependence of A,/(e; ip^) on e within the region 
|e| < ojd near iJj? is supposed and we assume that the pairing is mainly of s-type, 
that is the harmonics with J — Q give dominant contribution into ( 5.J ) . Substituting 
( 4.39 ) into ( ^ ) and recalhng ( 4.32 - 4.36 ) we then have 

;iPs,iPm) 



A(ip,)|i + ^T, 4!i(^*;i 



Tc ^ |i[)|^(e*;ips,i?'m)^'q(iPm)A(ip„)+ 

-L[,'^(e*; ip„ ip,„)Wq(ip,)A(ip,)} , (5.9) 



and this equation does not contain anomalous contribution due to disorder scatter- 
ing explicitly. 

If we introduce an auxiliary quantity S(ips) by the expression 

S(ips) = - Z{ips)\T(3 + x{t^Ps)t3 + Z{ips)A{ips)Ti, (5.10) 

we then can write for the isotropic superconductor the following set of linearized 
Eliashberg-type equations 

ips[l- Z{ips)] = -Tc ^ L[,^^(ips,ip™), (5.11) 
q,ip™ 

(x) 



Xiips) = -T, ^ 4^4(ip,,ip™), (5.12) 



q,ip™ 

^(iPs)A(ips) = Tc ^ Lj)^^(e*;ips,ip„)^'q(ip„i)A(ip„)+ 

qap™ 

+r, ^ 4'^(e*;ip„ip™)^q(ip,)A(ip3). (5.13) 
q,ip„^ 

Being considered as equations on auxiliary self-energy S(ips) these equations may by 
depicted in a form formally coinciding with that shown in Figs. 5 and 6, except that 
picking out of isotropic contributions is meant implicitly. We see that the anomalous 
disorder-scattering contribution to the self-energy can be eliminated owing to the 
Ward's identity for the disorder part of the self-energy similarly to the case studied 
in the previous Section and the remaining explicite dependence on disorder is due to 
renormalized EPI vertices involved in the equation for auxiliary self-energy E(ij5s). 
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The fact of the Ward's cancellations of anomalous contributions coming from 
disorder scattering and disorder self-energy contributions to the gap function man- 
ifests itself also in the so-called Abrikosov's identity.c3 To establish this identity, we 
consider the linearized <f>^{ips) contribution, 



i'lS^Ps) - ^{/j^i^(ips)*q(ips)^q(ips)7q(ip.)Aq(ip.), 



(5.14) 



to the function (j)-^{ips). Bearing in mind isotropic s-type pairing and making use 
of the expansion over harmonics ( 4.3C ) one obtains 



i>U^Ps) = C^iij(ip.)*q(ip.)^q(ip.)7q(ip.)A(ip.), (5.15) 



where we have retained a residual k-dependence of the irreducible four-leg vertex. 



Recalling (4.161) we have 



A(ip.) 



7k(iP«) 



(5.16) 



Substituting in this expression the definition (4.2) for 7j^(ips) and using the time 
reversal symmetry we arrive at the Abrikosov's identity 



</'|,(iPs 
A(i^ 



ilml]j,(ips 

ips 



(5.17) 



in the form proposed by Lustfield, so that this identity is not specific for BCS-like 
models but holds also for the case where pairing is due to EPI and irrespectively of 
the strength and type of quenched substitutional disorder. 



6. Qualitative Analysis of the Reduced Isotropic Equations 

General formal equations ( p.33 ,2.41- 2.43| ) from which one may obtain the explicit 
expression for the superconductive transition temperature Tc, are too cumbersome 
to work with directly. However, if certain approximations are made when calculat- 
ing isotropic contributions to S(ips), qualitative analysis of the reduced equations 
(5.11-5.13) becomes possible. Confining ourselves hereafter by the binary alloy 



model (2_^) introduced in Section 2 we present here general formulas for Tc in the 
case of weak and intermediate EPI coupling and for the renormalization of the EPI 
vertices by disorder scattering. 

To pick out isotropic contributions to T,{ips), let us consider terms of order 
M of the self-consistent perturbational series in V contributing to S(ips). The 
isotropic part resulting from a particular M'th-order diagram can be represented 
in terms of the functions G,j{e;ip,s) and Ajj' (ee'; iwn) which are the coefficients of 
the expansions over Fermi surface harmonics of the one-particle Green's function 



Gj^(ips) and of the function A]jq(ia-'„) (defined by (4.11)) respectively. 
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There arise also combinations of the form 

^o"'W E Fj,{k)Fj,{p)FjMFjA^-P + ci) 
kp.q 

6{e - t^)6{e, - tp)S{e2 - tq)S{es ~ p^^) (6.1) 

and the like containing more complex combinations of (5-functions and Fermi surface 
harmonics i^j(k). Quite convincing check of the statement made is to write down 
the expression for the isotropic part of S(ips) within simplest alloy approximation, 
that is VGA, but corresponding details are mainly of auxiliary character and we do 
not present them here. 

Our further approximations are formulated as follows. 

Firstly, we suppose that approximation adopted for Y.^(ips) is such as not to 
disturb drastically the k-dependence of G-^^{ips) comparatively with that of the 
bare one-particle Green's function, at least for the narrow energy region |e| < lud 
near Ep, what implies that S]jj(ips) is almost k-independent in this region. This 
holds exactly for all energies within arbitrary single-site approximation to S^(ips) 

and approximately for E^~^'*(ips), to a good extent however. More exactly, we 
may suppose that S^(ips) « E(f]^;ips), with ij^ being "bare" electron spectrum, 
and in particular we may have I]^(ips) — '^{^Ps) within an arbitrary single-site 
approximation. It is also assumed that dominant contribution to Ajj' (ee'; iaj„) 
results from harmonics with J = J' — 0. Hence, we have the following substitutions: 

Gj{e;ips) Go{e;ips)Sjo, (6.2) 

Xjj,{ee';iuj.m) Aoo(ee'; iwm)(5jo^.7'o- (6.3) 

In this situation combination (|6.lD simplifies 



N^\e) ^ S{e~t^)6{ei-tp)S{e2-tcOS{es-t-i^_p^^), (6.4) 
k,p,q 

and so do the like. 

Secondly, the dependence of Aoo(ee'; iw„) on e and e' is supposed to be sufficiently 
weak in the region |e|, |e'| < ujo and therefore may be completely neglected. This 
leads to the replacement 

Aoo(ee';i^m) — > Aoo(e*e*; it^m) = A(iu;m), (6.5) 

where e* is some point within the energy interval specified. 

Thirdly, we decouple expressions of the type (|6.4|) and the like in a simple manner 



No\e) ^ <5(6~tk)<5(ei-tp),5(e2-tq)<5(e3-ik-p+q)^ 
k,p,q 

N^\e) E '5(e"tk)<5(ei~tp)<5(e2-iq) JE^(^3-<k+q-p' 
k,p.q [ P' 

= No{ei)No{e2)No{e3). (6.6) 
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In these approximations one-particle Green's function enters the analytical expres- 
sions for diagrams in the combination 



(6.7) 



where the use was made of ( 4.30| -4.36). What happens when accepting these sim- 
plifications is that the vertices become purely local quantities. They contain only 
site-diagonal contributions, so that these approximations may be justified for the 
case where site-off-diagonal contributions carrying all remaining k-dependence of 
the vertices are not important. The latter holds true, in particular, when the lo- 
calization of electrons by disorder is absent or at least weak for the energies within 
the energy interval of interest. 

In any respect, these simplifications are fully transparent for single-site approxi- 
mations to the alloy self-energy where all irreducible vertices are purely local quan- 
tities and as for more complex non-singlc-site approximations we suppose the va- 
lidity of them, of course with the remarks above in mind. Note also, that the 
other way of introducing the approximations proposed is to use the limit of infinite 
spartial dimensionality d — oo^@'0 which has become very popular last time 
when discussing problems of the theory of strongly correlated systems .Ell^'Eiliio 
The problem of the renormalizaition of EPI by disorder may then be formulated 



in terms of the Holstein model^ in the limit d — cx£3 and with the disor der con- 
tributions added. This problem apeares to be a single-site problemESHcJ and the 
k-dependence of both the self-energy and the vertex corrections is absent from the 
very begining. 



And at last fourthly, the expression (6.7) can be represented as 



G{ips) = - V /de4^^{ip,ro + eT3 + A(ip,)Ti} 



= - de 



pi + 



{ipsTo + era -I- A(ips)ri} 



(6.8) 



where we have expanded Giips) up to the first order in A(ips) bearing in mind 
equation for Tc and introduced renormalized one-particle state density 



k 



(6.9) 



and A-^{e) is defined by (4.18 4.19). Choosing then Ep — and assuming weak 
dependence of N{e) on e in the interval |e| < lod we can write 



Giips) = -i7riV(0)Sign(p3)ro - ^^7r7V(0)ri. 

\Ps\ 



(6.10) 



It is this simple frequency dependence of normal part of G(ips) which enables, along 
with approximations (|6.2-3.6), to analyze qualitatively the equations ( |3.1l[ - pTT^ ). 
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Key idea of such analysis is this. After using approximations ( |6.2| - |6.6D and ( 6.10| ), 



decouphng of dependence on frequency and disorder potential V and concentration c 
proves to be possible order by order in the perturbational expansion. The equations 
can then be represented as 

^-im/"/3(^'^^(0)y)7r7V(0)A(ip, - ip^)Ca(ips)C/3(iPm), (6.11) 

where ip^-dependence of Ca(jPs) is known and Ca(}Ps) possesses one of the following 
simple forms 

1 c- ^ A Sign(ps) A(ip^) 

1, Sign(ps), — ■ — ■ — , — — — , (6.12) 

\Ps\ iPsl 

and the factors /a/3(c, 7r7V(0)y) in ( |6.1l[ ) are wholly determined by the approxima- 
tion used to calculate E'^(ips), and a, /3 — 1, . . . ,A. 

We consider this possibility in detail for Fock contributions to the self-energy 
first. 

Characterize diagrams of order M by three numbers N, L, Q, with N, L and Q 
being the numbers of successive disorder scattering processes before, between and 
after two successive EPI processes respectively (as we shall see shortly only a parity 
of N, L and Q is of importance). Obviously, the number of full electron lines in the 
A/'th-order diagram equals M + 1, the number of T3 Pauli matrices is M -|- 2 and 
N + L + Q = M. We also refer to the subset of graphs with triple NLQ fixed as to 
iViQ-family. 

T — I I • I — I 1 — I • J ; — 7 

N \ \ \ \ II ; / 
\ \ \ \ II / / 

NLQ 



Fig. 9. General structure of the Fock-type diagram. Scheme of AfLQ-partitioning. Broken line 
represents a single-site scattering po 
other elements see previous Figures. 



represents a single-site scattering potential; the oval, the quantity £^J^'q(c). For meaning of the 



Within the approximations (6.2^.6) and (6.10) the contribution of M'th-order 



self-energy diagram reads (see Fig. 9) 

B^nIqV'T, E ^('P^ - 'PrnH {G{}Ps)T:,f {G(ip„.)T3f + ' iG{ips)T3f (6.13) 

where B^^^q{c) is a product of the renormalized cumulants for the diagram {D) 
given. 
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To proceed with the calculation of normal part of S(ips) we note that diagrams 
with M being even or odd contribute to Z{vps) or x(iPs) respectively. Indeed, graphs 
with even (odd) A/ contain even (odd) number of T3-matrices, normal contribution 
to G(\ps) is proportional to unity ro-matrix (see expression ( |6.10 )), therefore the 
contribution of such graphs to S(ips) is proportional to tq (r3)-matrix. 

In turn, we further classify graphs contributing to Z{\ps) with the parity of N , 
L and Q. There are only four combinations for the decomposition of M even on iV, 
L and Q with definite parities: 

eee, oeo, eoo, ooe, (6.14) 

and e (o) means even (odd). 

Consider for example eee-type contributions. After inserting the normal part of 
G{ips), which has the form — i7riV(0)Sign(ps), general expression ( |6.13D reduces to 

5^%(c) (-i^)""' E - ip„0T-3'''+'Sign^+'^(Ps)Sign^+i(P™), (6.15) 
ip, 

where 

V^t:N{Q)V, \{iujn) = N{Q)\{iuj^). (6.16) 
Then, for eee-type contribution one has 

- ^T,Y,\(ips - ip„)Sign(p,„). (6.17) 

Analogously, for other parity combinations, we have 

-iB^%(c) (-iV^*^) ^T, k^Ps - iPm)Sign(p„), (oeo) (6.18) 

-i5wLQ(c) {-iV^') ttT, M^Ps - ip,n)Sign(p„), (eoo, ooe) (6.19) 

and, as we see, the precise form of frequency dependence for the expressions under 
sums is completely determined by the parity of the numbers N, L and Q. 

Putting everything together over all orders in V we then have the contribution 
of the Fock diargams to the equation on Z{ips) of the following form 

(c, V) + foeo{c,V)]7:T,Y, 

^{^Ps ~ iPm)Sign(ipm) — 

- ipm)Sign(ips), (6.20) 

and the coefficient functions /a/37 i^^ ^) 

U,ic,V)^ Y Y.BNLQ{c)[-iv) . (6.21) 

NLQ={al3-f} {D} 
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Here B^j^^q{c) is a product of cumulants for a particular graph (-D), inner sum runs 
over a set {D} of the Fock-type graphs with triple NLQ fixed, that is over the 
A^iQ-family, and outer sum is over all A^LQ-families with definite parity {a(3j} of 
NLQ triple. 

Similarly one obtains the Fock contribution to xiips): 

- iPm)- 

-i |/eeo(c, V") + /oee(c, t^)| Sign{ips)nTc^ \{ips - ipm)Sign(ip„). (6.22) 

Consider the anomalous contributions to I](ips). A particular M'th-order graph 
for the normal part of Y,{ips) generates M + 1 graphs contributing by successively 



replacing each normal state electron line with anomalous part of (6.10) which is 
proportional to ri-matrix. These graphs may be further partitioned into two com- 
plementary sets. The first set contains graphs where anomalous part of G{ips) is 
between two successive EPI processes; the second, all the graphs remained. 
For the first set we have 

B^^lgi^) (-i^) M^Ps - ip™)Sign^+«(p,) X 



,,n- L/ -.^(iPm) I N+S-1 L+Q+l-S ( no\ 

X Sign (p™)^—p <j 2^X3 ^ riT3 ^ !■ , (6.23) 



(L + l 



.5=1 



and the expression of the form 

> \ M 



A(ips), 



^re^A(ip.-ip„)Sign^+«-i(p,)Sign^+i(p,„)- ,, 

( N Q >| 

ls=i S=l J 

for the second set. 

Let now M be odd. Pauli matrices enter such diagrams in combinations 

(T3)°ri(r3)^=ir2, (6.25) 

(t3)Vi (t3)° - -ir2 (6.26) 

(here e (o) means even (odd) power respectivelly) , leading formally to non-zero T2 
contributions to anomalous part of Yiiips). However full contribution of such graphs 
vanishes identically owing to usual gauge symmetry under phase transformations. 
This may be established on rather general grounds for the equations ( p. 30 -2. 32 , 



2.41- ^.43 ) without any computational simplifications accepted when evaluating con- 



tribution of a particular graph provided the approximation to S^(ips) used possesses 
a property of being conserving. But the last is precisely what we deal with. 
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Now we check that withhi computational shiiphfications (|6.2-3.(: ,6.10) vanishing 
of r2-contributions is preserved. Consider a particular graph of asymmetric NLQ- 
family, N ^ Q. There is always a graph in the QLiV- family which is symmetric to 
initially given one by reflection with respect to a vertical line since full contribution 
to S(ips) possesses this symmetry. Both graphs generate M + 1 contributions to the 
anomalous part but with reversed order of Pauli matrices relativelly to each other. 
Therefore whole conribution of two asymmetric NLQ- and QLN- families vanishes 
identically (recall Eqs. ( 6.25 , |6.26|) ). Analogously, this vanishing takes place for the 
symmetric iVLQ-family, N — Q, but within this family, because graphs with the 
symmetry specified belong now to the same family: anomalous graphs of the first 
set cancel within the set and so do graphs of the second set. Thus we conclude that 
such contributions vanish identically order by order in the pertubational series. 

For even M, after partitioning on NLQ we have for the anomalous part of S(ips) 
the contribution which for the first set reads 

-<2q(c) (-i^)'' vrT,^ A(ip, - iP™)^^, (oeo) (6.28) 

and for the second set 

B^^lgic) (-i^)'' tT,^ A(ip, - ip„)^^Msign(p,)Sign(p„) (eee) (6.29) 



B'}^lQic){-^Vy' 7TT,J2H^Ps-^Prn)^l^- (eoo, ooe) (6.30) 



All other pariry combinations from (S.14) give zero contributions to the anomalous 
self- energy. 

Again, summing anomalous contributions over all orders of perturbational ex- 
pansion we obtain the contribution of the Fock diagrams to the anomalows part of 
the self-energy S(ips): 



{/eee(c,y)-/oeo(c, V)] TtT^ Y, 



A(ip„) 



[feooic, V) + fooe{c, V)^ ^py^TrTc ^ ~ iPn 



+2feeeic,V)^I^^TTTcy^ X{ips - ipm)Sign(ips)Sign(ipm), (6.31) 

\Ps\ ^ 

1Pt7 



where functions jap-^ are defined by Eq.(6.21). 

Now we briefiy consider evaluation of the Hartree-type contributions to the self- 
energy E(ips). The evaluation of these contributions follow closely same roots as 
above, and minor modifications are only necessary. 
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Fig. 10. General structure of the Hartree-type diagram. Scheme of AfLQ-partitioning. 



General Hartree-type diagram is shown in Fig. 10. 

We characterize Hartree-type diagrams by three numbers N, L and Q, where N 
and Q have the same meaning as before and L is the number of successive disorder 
scattering processes on a fermionic loop. Analitical expression for Hartree-type 
contributions has the form 

- 5^%l^*'r,^A(0)T3 (G(ip.)r3)^Sp{(G(ip™)T3)^+^} {Giips)T,f , (6.32) 

^n'lq i'^) product of the cumulants for particular diagram (D) with the triple 

NLQ fixed, and diagram (D) belongs now to the set {D} of Hartree-type diagrams. 
Then making the approximation ( ^.IC ) we have the contributions to the equations 



on Z{ips) and xi^Ps) of the form 

2i {/eoo(c, y) + /ooe(c, 1/)}A(0)ZiSign(ip,), (6.33) 

2i {/eoe(c, V) + foooic, V)} \{0)U (6.34) 



and to the anomalous part of the self-energy: 

U = ttT^ 1, (6.36) 

\Pb\<W 

and we have introduced the cut-off W, with W being of order of the bandwidth. 
Note also that the functions fap-y in the expressions ( 6.33 - 6^ ) are defined by 



- 2 {f,oo{c, V) + fooeic, V)} ~X{0)U^I^. (6.35) 

Here 



Eg. ( 3.21 ) where inner sum runs over Hartree set of diagrams. However the results 



of the summation for both Hartree- and Fock-type diagrams appear to be the same 
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as may be demonstrated by strightforward manipulations using the method of gen- 
erating functions. Moreover, it is not even necessary to calculate the coefficient 
functions (7.21) for Hartree-type diagrams because Hartree-type diagrams do not 
explicitly contribute to the final Tc expression within the approximations used in 
this Section. 

Finally, taking into account ( p^ , |6.22| , |Oll ) and ( |6.33| - |63^ ) one obtains the 
following set of the equations: 

ip,[l - ZCips)] = 2i [feooic, V) + fooe{c, V)^ A(0)Z^Sign(ip,) 

-i {/eee(c, V) + foeo{c, V)] TtT, ^ A(ips — ipm)Sign(ipm) 

-i [feooic, V) + fooeic, F)} TtT, ^ \{ips - ip™)Sign(ips) (6.37) 

Xiips)] = 2i [feoeic, V) + foooic, V)} ~X{0)U 
-i [foooic, V) + feoeic, V)} irTe ^ A(ip, 

-i{/eeo(c, V) + foee{c,V)]Si gn{ips)TrTc ^ A(ips — ipm)Sign(ip,„) (6.38) 



A(ips 



^ips)A{ip,)^-2[feooic,V) + fooeic,V)}~XiO)U 
[feeeic, V) - feee{c, V) | TrT^ ^ Xiips - iPm) ^^^'^"'^ 



\Pm\ 

[feoo{c,V)+fooe{c,V)\ ^^^^T, ^ A (ip, - ip„0 + 

2foeo{c, V) ^j-^^^^ ttTc ~X{ips - ip,„)Sign(ips)Sign(ip„i) (6.39) 

\Ps\ ^ 

Having been solved, the Eq.( |6.37 ) yields for Zivp^ the expression 

Z(ip,) ^ l-2{/eoo(c,F)+/ooe(c,F)}A(0)ZY^ 
+ {/eee(c, V) + /oeo(c, V)\ ^TtT, ^ A(ip, - ipm)Sign(ip^)Sign(ij>s) 

' ' i'Pm. 

+ [feoo{c, V) + /ooe(c, V)\ ^TtT, ^ A(ip, - ip™). (6.40) 

Substitution of ( 6.40 ) into ( 6.39|) then leads to the equation for the gap function: 
A(ip,) = vrT, ^ A^//(ip, - ip™)^^ 
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-^l^7TT,J2>'""i^Ps - ipm)Sign(ip,)Sign(ip„), (6.41) 



where 

A^^-^(iw„) = k(c, f )7V(0)A(ia;„) 

= K(c,f)A(ic^„), (6.42) 

k{c, V) = feeeic, V) - feoe{c, V), (6.43) 

with X^-^-^{iuJn) being EPI kernel renormahzed by disorder scattering and k{c,V) 
being the renormahzation factor depending on the approximation used for disorder 
self-energy. As we have stated earlier Hartree-type contributions do not explicitly 
appear in the equation for the gap function (3.41), which is the result of compen- 
sation between Hartree contributions to the function Z{ips) and those of to the 
anomalous self-energy (/>(ips). Of course, such compensation does not take place in 
generalli§Olla, and is rather due to the character of the computational simplifica- 
tions used, of which the local approximation (e.g., the limit d = oo) is principal. 
Comparing ( 6.41 ) with the standard result for A(ips) which rcadta 

A{ips) = ttTc A(ips - iPm)-^r-^T^- 

^ \Pni\ 

^P^ttTc X{ips - ipm)Sign(ips)Sign(ip„j) (6.44) 

\Ps\ ^ 

we see that the two equations differ in the only respect: the first contains A'^-^-' (ict;„) 
and the second, A(ictj„); the form of the equation remains intact. Then using square- 
well approximationP for Aoo(ips ~ ipm) in the first term of Eq.(6.41) 

Aoo(ips - ipm) — > XOiuji - \ps\)0{uji - \pm\) (6.45) 

and a slightly different version in the contribution which comes from the expression 
for Z{ips) 

Aoo(ip6 - ipm) — > X9{uji - \ps -Pin\) (6.46) 
we immediately obtain for 

r 1 + X^^f 1 

Te = l.liuD exp <^ \ (6.47) 



X^ff = k{c, V)N(Q)X,, a, = 2 / df7^^^^°^^^ (6.48) 



where 

/•oo 2 p 

\<^ff = K.(r V\N(n\\u \,. = 9 I HQ 

n 

The Eq. (|6.4l|) enables one to consider more general situation where, along with 
EPI mechanism, another pairing mechanism exists which will be called, for definite- 
ness, "excitonic". Here we treat only the case of weak coupling in both EPI and 
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"excitonic" channels. The corresponding interaction may originate from electron- 
electron interaction, but shoul d n ot constitute its main part. For instance, one 
may think of the Weber modclEll where strong Mott-Hubbard correlations in the 
d-system and pairing via crystal-field excitations are considered. 

In this situation "excitonic" channel leads to the appearance of the second scale 
in the effective interaction kernel and to the logarithmic corrections to \'^^^ . Within 
square- well model for the interaction kernel, Aoo(ips — ipm) may be approximated 
as 

Aoo(iPs - iPm) > \l9{uJi - \ps\)9{uJi - + \29{UJ2 - \Ps\)0{l^2 - bml), (6.49) 

where loi and L02 {^i < ^2) are characteristic phonon and exciton energy scales, Ai 
and A2 are corresponding coupling parameters of dimension energy. Substituting 
( |6l9| ) into ( pill ) we obtain 



= 1.13cjiexp|-^| (6.50) 
with A^-^-'^ defined by an expression of the form 

^ Ai + ^ (6.51) 

1 - A2 In |a;2/t^i| 

and 

\i = k{c,V)N{0)\i, A2 = K(c,y)7V(0)A2. (6.52) 



The expression (3.51) retains its usual form, but up to a renormalization of bare 
couplings by disorder scattering, which is absorbed in the factor k(c, V). 

Let us consider now disorder renormalizations of dilute paramagnetic impurity 



contributions to the gap function and Tc equations. In this case the equation (6.41) 
has to be changed to 

A(ip.) = ttT,^ A-//(ip, - ip„)^^- 
^^ttT, V X-lfiSps - ip„OSign(ip,)Sign(ip™) - 2fjf^^, (6.53) 

where 

7P = cpV];S{S + 1), TP = iV(0)7P, (6.54) 

t'^// = K(c,y)7>, (6.55) 

with cp being the concentration of the paramagnetic impurities, Vp being param- 
agnetic impurity scattering potential, S being the impurity spin, and some details 
of the derivation of ( |6.53 ) are given in Appendix B. We see that, analogously to the 



situations without paramagnetic impurities, effects of non-magnetic disorder lead 
to the replacement of bare couplings by the renormalized ones in the gap function 
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equation and the renormalization factor is the same for all bare couplings. The gap 
function equation preserves its form (compare with the case where non-magnetic 
disorder is absent), and the results for the critical concentration of the paramagnetic 
impurities and for the Tc suppression may be obtained by standard means.Q Thus 
we have the equation determining the critical concentration ric of the paramagnetic 
impurities 

7c = (1 + X'ff)T,o/2.26, (6.56) 

where Tco is the critical temperature with •j'^f^ = 0, jc = ncVpS{S + 1), and an 
approximate expression 

\u('^)=^{a + l/2)-^{l/2), (6.57) 



gives the Tc suppression by dilute paramagnetic impurities, 



eff eff rp 

=o.l4^±^ (6.58) 



7rT,(l + A-//) ■ 7c T, 
is pair-breaking parameter and ^'(2), the digamma function. 

7. Calculations of the Renormalization Factor within Conserving Single- 
site Approximations to the Alloy Self-energy. 

In this section we demonstrate that within an arbitrary conserving single-site ap- 
proximation determining the self-energy E'^(ips) in the normal phase of a super- 



conductor the renormalization factor k{c,V) (6.43) can be expressed in terms of 



suitably defined generating function for S'^(ips). To this end we recall that the co- 



efficient -B^^q(c) in ( |6.21 ) is nothing but a product of all the cumulants resulting 



from disorder-vertex contributions for the diagram (D) given, and these cumulants 
are the same which enter I]'^(ips), for the EPI renormalizations are consistent with 
the approximation to S'^(ips) by construction. 

Let us introduce Q-function, generating function connected with the number of 
successive disorder scattering processes, by the following expression 

Qiz^,z„z,) = J2 Y.^NlQi^>i4^l (7.1) 

NLQ {D} 

Here B^^q(c) is the product of the cumulants for the Fock-type graph (D) having 
the triple NLQ fixed, inner sum runs over all graphs of A^LQ-family and outer sum 
over all positive integers NLQ. Clearly, the functions fai3-y{c,V) can be expressed 
through Q{zi, Z2, Z3) by picking out even and odd parts of the latter with respect 
to each of its three arguments, that is 

faf3.yizi, Z2, Z3) ^ ^ ^ VlV2V3QiVlZl,V2Z2,mZ3), (7.2) 
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where rji — ±1, (i — 1, 2, 3), and 

/a/37(c, V") = /a/37 (2:1, 22,23) I 



(7.3) 



According to its definition the function (3(21,22,23) is the product of analogous 
functions for vertex corrections 



Q(2;i,2;2,23) 



j (21, 22)^3(22, 23) {1 + (56(21, 22, 23)} , 



(7.4) 



where we have introduced (J-functions for irreducible three- and six-leg vertices. 
(^3(21, 22) can be further expressed through (54(21,22) function, generating Q- 
function for four-leg vertex, as 



^3(21, 22) = 1 + (53(21, 22)(54(2l, 22), 
(53(21,22) = (1 - (54(21,22))"^ . 



(7.5) 
(7.6) 



Now we find the relations between the disorder self-energy E^(ips) and gen- 
erating functions introdiced above. From Eqs.(C.4, C.lC , |C.12 ) of Appendix C we 
have 



C/[G,G'] 
C/[G, G', G"] 



VP2{VG), 

V^Pi{VG,VG'), 

V^Pq{VG,VG\VG"), 



where 



P2(2) = ^q„(c)z", 



(7.7) 
(7.8) 
(7.9) 

(7.10) 



n=0 



and g„ (c) are renormalized cumulants for the single-site approximation being used. 
(Note that we reffere to the generating functions introduced in Appendix C as to 
generating P-functions.) 

Generating functions for four- and six-leg vertices can be expressed as 



^4(21, 22) 



Az- 



2ti\Jc (z-zi)(z-Z2) 

P2(2l)-P2(22) 



21 - 22 



P6(2l,22,23) = ^ / d 



^2(2) 



{Z - Zl){z - Z2)(2 - Z3) ' 



(7.11) 
(7.12) 



where contours G in the integrals encircle only the peculiarities of the denominator 
of integrands. Noting that the irreducible vertex of 2n'th order has n potential lines 
more than full electron lines and using (7.8,7.9) one obtains the relations: 



(54(21,22) = ZiZ2P4(2l,22), 
?6(2l ,22,23) — ZiZ223P6(zi, Z2, 23). 



(7.13) 
(7.14) 
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Finally, we write down the expression for k{c, V) once more for reader convinience: 

k(c, V) = /eee(c, V) ~ foeo{c, V) . (7.15) 

So, we see that the knowledge of P-generating function for I]'^[G] is sufficient to 
calculate renormalization factor k{c^V). 

The expression (7.15) for k{c^V) is not quite convenient to work with directly 
and may be transformed further. To this end, we introduce, in tight analogy with 
the Bethe-Salpeter equation (2.24), the generating function Ti{zi, Z2) by 

Vi[zi,Z2) = Pi{zi,Z2) + ZiZ2Pi{zi,Z2)Ti{zi,Z2). (7.16) 

(7.17) 
(7.18) 



Then from (|7J,|7jj7JJ) and ( |7.16| ), it follows 

Q3.{Z1,Z2) = 1 + ZiZ2Ti{zi,Z2), 

and 



i{zi,Z2) 



Ti{zi,Z2) 



Pa{zi,Z2)' 

Using ( 7.14|7.17 ) we can rewrite expression ( [7.4| ) for Q{zi,Z2, Z3) as 

Q(zi, 22,^3) = + Z1Z2T i{zi, Z2) + Z2Z:ir i{z2, z:i) + ZiZ2Z3SQ{zi, Z2, Z3), (7.19) 
where 

6Q{zi,Z2,Z3) = Z2T4{zi,Z2)T4{z2,Z3) + Qaizi, Z2)Peizi, Z2, Z3)Q3{z2, Z3). (7.20) 

To proceed further we use the Ward-type identity which relates P-functions for ver- 
tices of different order. From (C.17 ,C.19) of Appendix C the function ^6(2^1, Z2, Z3) 
can be expressed as 



P6{zi,Z2,Z3) 



Pi{zi,Z2) - Pi{zz,Z2) 



Zl - Z3 



In its turn, from (7.16) one has 



^^4 (2:1, 22) = {Ti^{zi,Z2) + Z1Z2) 



(7.21) 



(7.22) 



and using then (7.16 7.22) one can obtain 



Pi{zi,Z2) -^4(^3,22) 



P4:{zi,Z2)Pi{z3,Z2) , 



T4{zi, Z2)Ta{z3, Z2) 
X {dT{zi,Z2,Z3) - Z2T4{Z1,Z2)T4{Z3,Z2)) {Zl - Z2), 



where 



dT{zi,Z2,Z3) = 



r4(zi,2;2) - Yi{z3,Z2) 



Zl - Z3 



(7.23) 
(7.24) 
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Substituting now the expression ( [7^ ) into ( [r2ll ) one has 

p , N P4(zi, 22)^4(2:3, 22) ^ 
-Pelzi, 2:2,23) = :p-7 7X 

(9r(zi, 22,23) - Z2Ti(zi, Z2)r4(z3, Z2)) . 



Recalling (7.18) and Eq.(7.20) for 5Q{zi, Z2t z^), we can write 



(7.25) 



(7.26) 



(5(5(zi,z2,z3) = 9r(zi, 22,2:3), 
so that Q{zi, 22, 23) becomes 

Q(Z1, 22, 23) = 1 + 2l22r4(2i, 22) + 2223r4(22, 23) + 2i 22239r(2l , 22 , 23). (7.27) 

As for the combination feee{zi, ^2, 23) — foeo{zi, ^2, 23), a little algebra with the use 
of (|7j) yields 

^ {Q(2l, 22, -23) + (3(-2l, -22, 23) + (9(21,-22,-23) + (5(-2l, 22, 23)} . (7.28) 

Noting the symmetry Q{zi, 22, 23) = (9(23, 22, 21) one has 



k(c,F) = ]^{Q{z,-z,-z) + Q{-z,z,z)} 



= -iV 



]{Q{z, ~z, -2) + g(2, 2, -2)} 
^{(9(21, 22, 22) + (9(21, 21, 22)} 

Substitution of ( |7.27D into ( |7.3lD leads 



(7.29) 
(7.30) 
(7.31) 



z\ = — zn =—\V 



k(c, y) = 1 + 2(zi-Z2)^ {2?(-Zl - 22)r4(2i, 2i) - 2fr4(22, 22) + 
+ 22122(21 - 22)^r4(2l,22)} = 

and on a line 21 = — 22 this expression further reduces to 

z\V^{z\,Z\) - 2zlz2T4{zi,Z2) + 2|r4(22, 22) , 



k{c,V) = 1 + 



{zi - Z2Y 



'zi = — Z2 =—iV ' 



(7.32) 



(7.33) 



k{c, V) = 1 + \^2ini^jm,,r4(2i, 22) \,^^,^=iy , (7.34) 

where Im^ means picking out "imaginary" part with respect to a complex variable 
2, which may be written as 

1 



lmj{. . . ,2, . . .) 



2i 



(/(■ 



•)) 



(7.35) 



and points stand for variables other than 2, and we retain explicit dependence on 
2i and 22 in (7.33,7.34) for further use when analyzing the renormalizations within 
CPA. 
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Now everything is ready to evaluate k(c, V) within certain widely known single 
site approximations. To be concrete we consider three such approximations: VGA, 
ATA and CPA. 

7.1. Virtual Crystal Approximation and EPI Renormalizations. 

For the case of VGA the self-energy reads0 

SycA(iPs) = cV + c(l - c)VGiip,) (7.36) 
and cumulants <?,i(c) are 

<?o(c)=c, gi(c) = c(l - c), g„(c)==0 n>l. (7.37) 
The generating function P2 (c) within this approximation has the form 

P2{z) = c + c{l~ c)z (7.38) 

and the use of ( fT.ll ) gives 

P4{zi,Z2)=c{l-c). (7.39) 



Then, from (7.16) and within the VGA accuracy, we have 

r4(zi, Z2) = P4{Z1, Z2) = c(l - c) (7.40) 



and from (7.32,7.34) it follows 

k(c,1/) = 1. (7.41) 

So, there is no disorder contributions to the renormalization of EPI coupling, except 
for usual ones through one particle state density. 

7.2. Renormalizations of EPI within Average T-matrix Approximation. 

For ATA approximation the expression for I]^j,^(ips) has the formEl 

^^^^(^^-)- l-(l-^W) - ^'-''^ 
General expression for (7„(c) within ATA reads 

<7„(c) = c(l-c)", n>l (7.43) 

and the generating function P2{c) is 

P2{z) = c[l - {I - c)z\-\ (7.44) 

Then, similarly to the previous case, we obtain 

P4[z^,Z2) = c(l - c) [1 - (1 - c)zi]-i [1 - (1 - c)z2r' (7.45) 
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and to required accuracy 

r4(zi,Z2) = P4(^l,Z2), (7.46) 

so that in this case the resuh for k(c, V) reads 

c(l - cYV'^ 

n{c, V) = l+ ^ ' \ . (7.47) 

(l + (l-c)2F2)2 

Let us consider now the so cahed spht band hmit, V oo. Supposing that the 
chemical potential lies in the lower A-subband, that is 

lim Nb{Q) = 0. (7.48) 

V->ao 

Introducing projected one particle state densities for A- and B-type atoms 

N{e) = [l-c)NA{e)+cNB{e), (7.49) 

we find in this limit that 

k(c,cx))= lim k(c,T>) = (1 - c)-i (7.50) 

V^oo 

and effective EPI coupling becomes 

X^ff = k{c,oo)N{Q)\ = \Na{Q). (7.51) 

Such a limiting behavior is fully transparent and might be expected on physical 
grounds. Indeed, suppose V is so large that electrons can hop and interact with 
phonons only on A-type atoms (the probability of finding an electron on B-type 
atom is exponentially small in this case, so that one may completely neglect this 
process). Then two effects arise. The first is the renormalization of the electron 
state density for A-type atoms by both B-type impurities and EPI. The second is 
the renormalization of chemical potential which is due to the fact that B-type atoms 
are no longer accessible for the electron motion. The electrons can now propagate in 
those states which belong to the renormalized band for A-type atoms and the pairing 
occurs for these states only. The temperature Tc depends solely on the properties 
of A-band states but not on the (macroscopical) size of superconductive cluster and 



this is the reason why X'^f^ obeys the limiting expression of the type (7.51) where 
the dependence on the concentration of B-type impurities and on the strength of 
disorder scattering enters implicitly through the renormalized electron state density 
Na{^) and, of course, throvigh the renormalized chemical potential. However the 
self-consistent version of ATA used here fails in discribing normal state properties 
of an alloy (see Ref.EH and references therein for critique on this approximation), so 
that more complex approximations have to be considered. 



7.3. Coherent Potential Approximation and EPI Renormalizations. 
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For CPA the equation for a self-energy E^p^(ips) readJ^l'Ei 

< t{ip,) 

where 



t(ips) = {Vrit 



^-'CPa) + i^Vi — ^CPA 



)Gii{ips)t{ips). 



(7.52) 
(7.53) 



We may introduce ^2(^) generating function by 

< t{z) 0, 

t{z)={l^-P2{z)) + {l^~P2(z))zt(z), 

where is a random variable with the distribution 

Prob{r; = 1} = c, Prob{?7 = 0} = 1 - c. 



(7.54) 
(7.55) 

(7.56) 



To calculate Pa{zi^ Z2) we use the Ward- type identity (|7.11| ). Then a kind of trans- 
formations used in Appendix A, Eq.(A.5), leads 



Pi{zi,Z2) 



< t{zi)t{z2) > 



I + Z1Z2 < t{zi)t{z2) >' 

and from Bethe-Salpeter type equation ( 7.16| ) one has 

^4{zi,Z2) =< t{zi)t{z2) > ■ 



Using (7.34) we can write k(c, V) in the form 

k{c, F) = 1 + f 2 < {lmt{iV)f > 



(7.57) 



(7.58) 



(7.59) 



A more familiar form arises if one uses general expression for k(c, V) in the form 
(7.33). To this end we rewrite the expression ( 7.33[ ) in the form 



~ _ < jzi + Zjtjzi) -Z2- Zlt{z2)f > 
^i^^^) - (zi - Z2)2 \z^=-.. = -iV , 

where the use was made of Eqs.( [7.53| ) and ( [7.57| ), and finally 

< [lin{z + z'^tiz))]'^ > 



k{c,V) 



\z = -iV 



(7.60) 



(7.61) 



(Imz)^ 

Within the approximation ( B.lOj ) one has 

ReG(+iO) = 0, ImG(+iO) = -i7riV(0), (7.62) 

with G{E + \Q) being analytical continuation of G{\ps) onto the upper half- plane of 
the complex energy plane, and the expression (7.61) becomes 

< ^'(0) > 



k{c,V) 



<iV„(0) >2 



(7.63) 
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where in the course of transformations we have introduced projected electron state 
densities as 

N^{E) = --ImF„(^), F^{E) = G{E) + G{E)mE + iQ)G{E), (7.64) 



with t{E + iO) being analytical continuation of t(ips), Eq.(7.53), onto the upper 



half-plane and G{E) being on-site one-particle Green's function. 



As one may see the formula (7.63) is exactly the same as in Refs.00 and obeys 



the expected limiting behavior Eq.(7.51) (the reasons for this are the same as in the 
case of ATA). 

8. Conclusions 

In this paper we have studied the electron system which includes both quenched sub- 
stitutional disorder and electron-phonon interaction. We have developed a rather 
general approach for constructing Eliashberg-type equations describing supercon- 
ductive transition in such systems and have analyzed a linearized version of the 
equations in a number of widely known "alloy" approximations to find explicit ex- 



pressions for T^. The examples are Eq.( 7.63 ) for CPA, which coincide with that 



given in Refs.uD, and also Eq.( 7.47 ) for ATA (along with general expression ( 6.47 ) 



for Tc and Eq.( 6.43| ) for the renormalization factor k(c, V) of course). 



Analogously, one may derive rather simple and compact expressions for in any 
other "alloy" approximation for normal phase and within more general models of 



quenched disorder including off-diagonal one. The structure of the equations (2.30- 



2.33 2.41- 2.4^^ ) remains intact irrespectively of the type of quenched substitutional 



disorder as well as the approximation for the "alloy" self-energy used (obviously, 
the latter must obey (/)-derivability criterion, Eq.( [3.44| )). 

We want also mark a general conclusion on extreme importance of EPI vertex 
renormalization by alloy disordering and on the Anderson's theorem for arbitrary 
type of quenched disorder treated within conserving approximations. As for results 
( |6^50| - |6^ ) and ( |6^55| - |6^ concerning the influence of disordering on logarithmic 



corrections to the effective EPI coupling and on dilute paramagnetic impurity sup- 
pression of Tc respectivelly, to our knowledge they have been presented in such a 
general form for the first time here. 

It would be interesting enough to generalize the approach presented to the case 
of non-singlet types of pairing, which seems to be quite possible. In this case, as 
one knows, the Anderson's theorem does not hold and the effects of the disorder 
influence on Tc are much stronger .i 

Let us discuss briefly the applicability of the results for Tc obtained to the "alloy 
analogy" for the Hubbard model. Seemingly, they are not valid literally because of 
very different structure of diagrams for quenched disorder and for effective disorder 
within static approximation for the Hubbard model, where disorder is not quenched, 
but annealed. In this case the "concentrations" of "alloy" components themselves 
are functionally dependent on one-particle Green's function and should be varied 
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too. At a diagrammatic level of discussion the difference among two types of the 
disorder is the appearance of graphs with loops contributing for the case of annealed 
disorder, while such graphs are absent for quenched disorder. But even if one 
neglects this difference (which may be quite reasonable, at least for not too strong 
on-site repulsion when distributions for fluctuating fields are weakly renormalized 
by electrons) there is still another reason for why the present results are not directly 
applicable to the "alloy analogy" for the Hubbard model. The matter is that in 
this case fluctuating fields contain both charge and spin components, so that the 
effective disorder is by no means nonmagnetic. 

One can also see that simple approximations of a type used by us earlieril are 
rather dangerous since a type of Ward's cancellations may take place just as in the 
case of Anderson's theorem and careless generalization of the present approach to 
the "alloy analogy" within functional integral treatments of the Hubbard model is 
required. 

Let u s strpRs once more the difference between our treatment and the approach 
by BelitzllaOlla Belitz, following ideas due to Anderson,uH used the superconduct- 
ing order parameter A^, cx cj,-^C^^ where i' and V are the states connected by the 
time reversal symmetry for a particular disorder configuration. At the same time 
we use k-dependent order parameter in a translationally invariant effective medium. 
Therefore both approaches are not equivalent from formal point of view: diagrams 
of the same form in ly- and k-representation need not to give the same contribu- 
tion and, hence, the same result for Tc- In particular, in the approach by Belitz 
the Hartree-type diagrams do not contribute to the anomalous self-energy unlike 
present approach. 

To establish the equivalence one has to solve very difficult, yet not investigated, 
problem whether the random quantity A^, is a self-averaging one. On the other 
hand, the order parameter used here seems to be directly related to the quasiparticle 
spectrum in the superconducting phase and consequently to observables. To our 
opinion, treatment of such important problems as, for instance, order parameter 
anisotropy in high Tc superconductors is more adequate within our approach. In 
any case, the order parameter used here is physically satisfactory, provided the 
quasimomentum remains well defined, that is, far from the localization transition. 
As for the essentially strong disorder (near the Anderson's transition), both the 
approaches apparently do not work and the choice of the order parameter in the 
intermediate region may be done empirically. 
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Appendix A The Ward identities 

In this Appendix we present the derivation of Waid's identities ( ^.41 ) and ( 4.15| ). 
The identities of this type was established in Ref.E3 within the CPA in connection 
with the calculation of transport properties of binary alloys. The relations between 
self-£nergy I]'^[G'], as given by_CPA, and irreducible four-leg vertex presented in 
RefEj was also studied in Ref.E^ with the use of diagrammatic technique. 

To begin with, we establish the connection between variation of S[G] caused by 
finite change in G. Such an identity was derived within certaiiL approximation and 
used when discussing the problem of localization by disorder .Ej 

Given arbitrary but finite variation in one-particle Green's function G 

G' = G + AG, AG — finite! (A.l) 

we consider 

< AT >=< -T'^ >, (A.2) 

where T'^ is given by 

T^iips) - {V' - ^'{ips)) + iV' - ^%ips)) Giips)T%ip,), (A.3) 

T'"^ is defined by similar expression, but with G and S replaced by G' and E' 
respectively, 

E' = S'' + AE, (A.4) 

and we impose the self-consistency condition < T >=< T' >= for E to acquire 
proper, albeit implicit, functional dependence on G. 
The following chain of equalities holds: 

=< AT" >=< T'^[G', E'] - T"[G, E] > 

= - < ((r'")-i - (r'^)-i) > 

= - < {{V" - E'")-i) - - E-^)-!) - (G' - G)) T" > 
= - < (1 + T"=G') AE(1 + GT") - T'^AGT" > 
=< VAGT" > - < (1 + T'"G')AE(1 + GT") > -AE. (A.5) 
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In explicit index notation the last line of Eq 



) reads 



[G',G]AG„ 



/3' 



(A.6) 



where, analogously to ( |3.2S| ), we have introduced full disorder four- leg vertex 

r^a;/3L/[G', G] =< T'^^^Tp^ > . (A. 7) 

The Bethe-Salpeter equation gives the irreducible four-leg vertex 

r^Q;/3,y[G', G] — C/^Q;/3,y[G', G] + U -^^^ i,[G' , G]G' ^' a' G f}' ly'T a' a;f3l3' [C , G]. (A. 8) 

Multiplying ( |A.6D by the combination G'?7[G',G]G and using the Bethe-Solpeter 
equation (|A.8|) we can write 



rQa';/3'/3[G', GJG'a'pGiy/J' AS^y — 
Uaa' -,0' I3[G' , G]G' a' fj.' G^' I3'T ^/ s-ju' [G', GJAG^^. 



(A.9) 



After substituting this expression into (A^) and using Eq.( |A.8[) once again one 
obtains the Ward's identity 

[G',G]AGa0. (A.IO) 
Putting now G' = G -|- AG, where AG tends to zero we immediately have 

^ = C/pa;/3.[G,G]. (A.ll) 



(A.12) 



By the definition (A. 7) the vertex function r[G', G] has the property 

^ ^ji,a;l3v[G' T G\ — r^jy.pQ [G, G'] 

which induces corresponding property for irreducible four-leg vertex 

Ufia-f3v[G' ,G] — Ui3u-iia[G,G']. (A. 13) 

Then the use of (A.ll) leads to the condition of (/)-derivability, and the Eqs. ( 3.26| ) 
which implicitly determine the self-energy S'^[G] can be represented in the form 
stated in Eq. (2.33) with some functional M^'^[G]. 

Note also that for any conserving approximation for I]'^[G] the condition of 
0-derivability holds automatically and the identity of the type ( A.IC ) may be es- 
tablished along the lines of Appendix of RefE3 and making use of the fact that 
complete families of self-energy terms contribute to E'^[G] in this case. 

Formal identity ( A.IO ) may be used to derive various relations between alloy self- 
energy Y.'^iips) and one particle Green's function Giips). For example, substituting 



S' = S^(ip3), S = S^(ip„ 



(A.14) 
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and 

G' = G'^iips), G = G,(ip™) (A.15) 



into ( A.IC ) we obtain 

E^J^{ips) - i;-^-^(ip,„) = ?7i;;„j(ips,iPm){G/„o-(ips) - Glnai^Pm)} (A. 16) 

which after Fourier transforming and using 

lmJ:i^{E)^^^[i:i^iE + i6)-i:i^iE-iS)}, (A.17) 

ImGk,(£;) = i{Gk^(i? + i5)-Gk,(£;-i5)}, 6^+0, (A.18) 
takes the form 

lmi:i^iE)=Y,U^^iE+,E~)lmGi^jE), E^=E±iO. (A.19) 

q 

This is nothing but an analogue of the unitarity condition for alloy type problems 
we deal with in the present paper. 

Appendix B Disorder-Induced Renormalization of Paramagnetic Impu- 
rity Tc Suppression 

In this Appendix we treat the problem of how Tc-suppression by dilute paramag- 
netic impurities is renormalized by non-magnetic disorder. Using the approach of 
Section 3 we derive the expression for the self-energy correction which is due to 
the scattering on the paramagnetic impurities, and then we apply the analysis of 
Section 5 and 6. 

Three assumptions are made. Firstly, paramagnetic impurities are distributed 
randomly with the concentration cp over the host and independently of the non- 
magnetic ones. Secondly, the disorder due to paramagnetic impurities is purely 
site-diagonal. Thirdly, both the concentration of paramagnetic impurities and the 
scattering intensity Vp on them are small, so that the relevant small parameter is 
cpVl 

The contribution to the Hamiltonian ( 2A_ ) due to paramagnetic impurities reads 
Hp = J2 {^rcJroC,; + 5+cJr_Q + Srcjr+C,} , (B.l) 

i 

where cj and C.^ are row and column spinors ( 2.10 ) on a site i, SI is z-component 
of the impurity spin operator of spin S", the matrices t+ and t_ are 1/2(ti + vr^) 
and 1/2 (ri — yti) respectively, and S~ have similar definitions, and 

Vt = ^pC, (B.2) 
Prob{G = 0} = (1 - cp), Prob{G = 1} = cp. 
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with being independent random variables. Up to the lowest order in cpVp, one 
has the correction to the effective potential V(ips) Eq.(|3.3[), 



SVapiiPs) = cVpS[S + l)gal3{iPs)5i^,i^, (B.3) 

and we have neglected uninteresting contribution of the first order in Vp which gives 
a shift of the chemical potential. 

The expression (B.3) can be written in the form analogous to the effective scat- 
tering potential caused by EPI 

5Vai3{\ps) = S^yKl^,ga'f3'{ips)K'^,0, (B.4) 

{K},.,^. = SuS.iTo, (B.5) 

S^Y = Sw = S-,YCpV^S{S +1) = Su,cpV^S{S + 1). (B.6) 

The procedure of Section 3 then yields the correction to the full self-energy T,[G, D], 
with the form 

SEf,ty{ips) = S^Y 

xG^,>a{iPs)Rla'i^Ps,'^Ps)Ga's{iPs)G^p{ips)R]jf^,{ips,iPs)G(}'^'{ips) }. (B.7) 

In this expression R'^p(ips,ips), is renormalized "paramagnetic impurity" vertex 
which obeys the equation 

^L'(iP^'iP"i) = ^la' + '^al3:,l3'a'{iPs:iPm)Gpf,{ips)Gf,'i3'{ip,n)K^f„, (B.8) 
-RL'fe,iPm) + UaP:,P'a'{iPs,iPm)Gi3f,{ips)Gf,>p'{ip,n)Rlf,,{iPsAPm)- (B.9) 



The expression (B^) may be depicted in a form similar to that for E'^^p'*[G, D] 
(see Fig. 5) with graphical representation of the combination cpVpS{S+l) replacing 
phonon lines and with the "paramagnetic impurity" three-leg vertices replacing EPI 
vertices. 

Further manipulations repeat literally the analysis of Section 5 and, after adopt- 
ing binary alloy model for non-magnetic disorder, we apply the qualitative analysis 
of Section 6 to the reduced equations obtained. 

Introducing 

7P = cpV^S{S + 1), 7P = iV(0)7P, (B.IO) 
one can write the expression for the correction to Z(ips) as 

SZ{ips) = TT^P {/eee(c, V) + foeo{c, V) + feoo{c, V) + fooe{c, '^^l ^ ^^'^^^ 

and for the contribution to xi^Ps) has the following form 

5x(SPs) = -i^7P [fooo{c. V) H- /eoe(c, V) + /eeo(c, V) + foee{c, V)] (B.12) 
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which has no energy dependence (within the approximations of Section 6, of course) 
and leads to an uninteresting shift of chemical potential only. 

The contribution of paramagnetic impurities to the anomalous part of the self- 
energy has the form 

= {feee{c,V) - feoo{c,V) - fooe{c,V)--3foeo{c,V)} (B.13) 

and the correction to the gap function becomes 

<5A(ip,) = -2j'0'^^j^. (B.14) 

\Ps\ 



In Eg. ( B.14 ) we have introduced 

Jp" = i^ic,V)jP, (B.15) 

where k{c,V) is defined by Eq.( |6.43| ). So we see that, up to the factor k{c,V), 
the contribution of the dilute paramagnetic impurities to the equations for the self- 
energy has the same form as in the case without alloy disordering. 

Full equations on Z(ips), x(ips) and A(ips) which take account of the paramag- 
netic impurity contributions read 

ip,[l - Z{ips)] = 2i {/eoo(c, V) + fooeic, V)} A(0)Z^Sign(ip,) 

X{ips - ip„i)Sign(ip,„) 

-i {/eoo(c, V) + fooeic, V)] nT, J2 H^Ps - ipm)Sign(ips) 

-i^7P {/eee(c, V) + ,Uoo{c, V) + fooeic, V) + foeo{c, V)} Sign(p,), (B.16) 

XiiPs)] = 2i {/eoe(c, V) + foooic, V)} A(0)Z^ 

-i |/eeo(c, V") + foee{c, V)^ Sign{ips)nTc ^ X{ips - ip„)Sign(ip„) 

{foooic, V) + /eoe(c, V) + feeo{c, V) + foee{c, V)} , (B.17) 



-ITTJp 



A(ips 



Z(ip,)A(ip,) = -2 {feooic, V) + fooeic, V)} A(0)Z^ 
+ {feeeic, V) - /eee(c, V)^ TtTc ^ Xiips - ip,n) ^^^^ 



I 

\Pr, 
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+ {/eoo(c, V) + fooeic, V)] ^^^T, ^ A(ip3 - ipm) 

~2/„eo(c, y)^^^r, > ^ A(ip, - ipm)bign(ip5jbign(ip 



-^7P {/eee(c, V) - feoo{c, V) - fooeic, V) - 3foeo{c, V)} . (B.18) 



Solving the equation (B.16) on Z{ips) and inserting into (B.18) one finds 

A(ipm) 



A(ip,)=7rT,^A^-^-^(ip, -ip„ 



\Pr, 



-^^nT, J2 X'^HiPs - ip„OSign(ip,)Sign(ip„) - 27^/^^^, (B.19) 

— the equation on the gap function, from which results foiithe critical concentration 
and for Tc~suppression follow by standard considerations 

Appendix C Generating Functions for Single-site Approximations 

In this Appendix we establish relations between suitably introduced generating 
function for arbitrary conserving single-site approximation for disorder self-energy 
T,'^ [G] and generating functions for irreducible disorder four-leg vertex and for higher 
vertices. The Ward's type identities for generating functions will be also derived. 

Generating functional for arbitrary single site approximation has the following 
analytical form 

W^G] = E E E In-lic)^-^^^^^^^ (C.l) 

and outer sum runs over lattice sites. For the corresponding self-energy one has the 
expression 

OQ 

S-^iG] = vY^lnic) {VGu{iPs)r , (C.2) 

ri=0 

where qn{c) is n-th order cumulant renormalized by the multiple occupancy cor- 
rections within approximation of interest. Particular choice of these cumulants 
determines wholly the approximation for the disorder self-energy. We introduce the 
generating function: 

00 

The expression for S'^[G] can then be written in the form 

Y.''[G]^VP2{VG). (C.4) 



Note also that generating functions we use here are connected with the number of 
Green's functions in the corresponding contributions to the self-energy £^^[0]. 
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Diagrams for the irreducible disorder four-leg vertex are generated by varying 
corresponding contributions to the self-energy and the expression for C/4 [01,62] 
reads 



00 CxD 00 



[/[Gi, G2] ^ v'j2 E E (vG,r ivG2r 

00 00 

= ^'EE9™+'+i(^)(^Gir(^^2r, (c.s) 



m=0 1=0 



where we have used the fact that for quenched disorder the self-energy I]^[G] does 
not contain fermionic loops and therefore the disorder four-leg vertex may be consid- 
ered as some functional of two independent functional variables. Generally speaking, 
we may think of disorder 2m-leg vertex as a functional of m independent functional 
variables by the same reason. 

Analogously to (C.S) we introduce the generating function 

CO 00 00 

n=2 m=0 1=0 

00 00 

= EE9™+'+i(c)2"2", (C.6) 

m=0 1=0 

and hence 

C/[Gi,G2] = V^Pi{VGi,VG2). (C.7) 

In general, for 2m-leg disorder vertex within arbitrary single-site approximation 
one can readly obtain 

00 00 00 m 

?7[Gi, . . . ,G„] = V^Y. E • • ■ E Sn^n,+-+n,^^+m J] (VG^P 

n—rn ni— n^— j — 1 

oo oo rn 

= ^"^ E • ■ ■ E 'Z^.+.-.+n^+m-l W n i^G.r , (C.S) 



m— n,„— j — 1 



so that 



cx) oo oo m 

^2m(^l, ■ ■ ■ , Zjn) = ^ qn-l{c) ^ • • • ^ 

n—m Til —0 — J — 1 

oo oo 711 

E'" E ^m+'-'+n^-Hm-lW JJzJ', (C.9) 



m— 71^,— j—1 



and then 



;7[Gi,...,G™] = y'"P2™(yGi,...,l/G™). (C.IO) 
We now express cumulant of the n-th order in terms of Cauchi integral 



58 A.O. Anokhin & M.I. Katsnelson 



The substitution of <7„(c) in the form ( p. 11 ) into ( p.9[ ) leads 

ni^O n^^O — ^ 

^ ^ ^ m I oo 

= — / AZP2{Z)\{ , 

27^1 Jc 2; - z. 



(C.12) 



where the contour of integration encircles singularities of the denominator only, and 
P2m(,zi, . . . ,Zm) posscss the property 

P2m{zi, ■ • ■ , Zra) — P2m{Z'Pl, ■ ■ ■ , Z-pm), (C.13) 

with V being arbitrary permutation of indices {1 . . . to}. 

The expression (C.12) relates generating function P2{z) with the corresponding 
generating functions for 2TO-leg vertices. In particular, we have 



and 



^4(21,^2) = J-T [ dz 
27ri Jc 

Pel^i, 22,23) 



P2{Z) 



dz- 



{z - zi){z - Z2) 
P2{Z) 



2tt\Jc {z - zi){z - Z2){z - Z3)' 
For TO > 2 rewriting the product in ( p. 12 ) as 



(C.14) 
(C.15) 



m ^ 



1 1 ™ 1 

— — 

Z - Zl Z - Z2 ^ \ Z- Zj 



Z2 — Z\ [Z — Z2 Z ~ Zi 



- ri -L J. r — ■ 



(C.16) 



J=3 



and substituting this expression into (C.12) we obtain the analogue of the Ward's 

3r 

(C.17) 



identity which connects the generating functions for vertices of different order 

P2{m-l){z2, Z3, ■ ■ • , Zm) ~ P2{m~l){zi: ^3, . . . , Zm) 



-P2m(2l, . . . , Z„i) — 

This equation gives 
and 



^'4(21, Z2) 



P6{zi,Z2,Z3) 

for m = 2 and to = 3 respectively. 



Z2 - Zi 

^ P2{Z2) - P2{zi) 
Z2 ~ Zi 

Pi{z2,Z3) - Pi{zi,Zz) 

Z2 — Zi 



(C.18) 
(C.19) 
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Yet another way to establish the relation ( C.l§| ) is to use the Ward's identity de- 
rived in Appendix A. Indeed, noting that within arbitrary single-site approximation 
irreducible four-leg vertex is a purely local quantity we rewrite (A. 7) as 

V{P2iVG2) ~ PiiyCx)} = VPi{VG2, VGi) {VG2 - VGi} , (C.20) 

where the use was made of ( |C.4|) and (C/?). Then substitutions 

VGi -> zi, VG2 ^ Z2 (C.21) 

immediately lead to the expression ( |C.18 ). 



